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® Fflr CAS- An finger, CAK-F-url^cSChriltcrte, 

Rcferenten, DMaklikcr 


Eberhard Lehmann 
LdvSoft, www , sna fu d & r - m i its 

Berlin August 2007 


Another new publication is: 

Gunter Scheu u.a., ABITUR - Priifungsaufgaben mit Losungen, Mathematik CAS, 
Gymnasium Baden-Wiirttemberg, STARK Verlag 

If you know about CAS- related publications then please inform me. Books dealing 
with the use of technology in math education in general are also interesting. Josef 
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Dear DUG Members, 

I am sorry for the delay delivering DNL#67. But I 
believe that its contents is worth some days over- 
time. We have again much more pages than the 
usual 40 - 44 pages. I didn't want to leave one of 
the articles for the next issue. Josef Lechner's con- 
tribution is one of the many ones which waited very 
long for publication, here it is at last. Nils Hahnfeld 
has been busy as ever and produced new packages 
for the handheld. As one of my final diploma pa- 
pers at university dealt with conics I was very in- 
terested in his CME ( Conics Made Easy). I could 
not resist making suggestions for including addi- 
tional features and a very intense communication 
has started. You can find a selection of this discus- 
sion together with a short presentation of CME. 
Differential Equations Made Easy will be presented 
in a future DNL. 

One highlight in this issue is doubtless Albert 
Rich's visionary paper. I am very happy and proud 
that Albert like so many of you keep on communi- 
cating within the DUG. It would be great if DUG 
members would follow Albert's invitation for sup- 
porting his vision. 

Another vision is to transfer as many Derive func- 
tionalities as possible to its intended successor, the 
TI-Nspire. The Tl-programmers promised to do 
their best and they are asking for the DERIVE- 
people's input. The DUG can play an important 
role in transfering the “DERIVE spirit 64 to the 
Nspire. (Maybe it is not a pure chance that Derive 
SPIRi t and NSPIRq have an important part of their 
names in common). Bernhard Kutzler presents a 
more detailed view on the Nspire. Next releases 
should contain one or the other of the DERIVE- 
features. Bernhard will continue in the next DNL 
with a special example to demonstrate how the 
various applications of Nspire can cooperate. 

At the occasion of a meeting with the TI- 
programmers some of us presented our most appre- 
ciated Derive features, which we would like to have 
implemented in future versions of Nspire David 
Sjostrand showed an excellent example of applying 
the ITERATES-command. This and other interest- 
ing findings in a triangle are collected in David's 
contribution in this DNL. He extended some ideas 
into 3D space (see next DNL). 


Another gem is promised for one of the next DNLs: 
A. Perotti joined the DUG recently. I remembered 
Alessandro's article on Groebner Bases in an ear- 
lier issue of the Inti DERIVE Journal. I wrote a 
mail and he immediately answered. He has adapted 
his DOS-Derive functions for DERIVE6 and wrote 
a full paper. He gave permission to publish it in the 
DNL and additionally promised to write an extra 
introduction. So we can hope that solving systems 
of nonlinear equations applying Buchberger's algo- 
rithm will not remain a „Black Box“ for some (or 
many?) of us. 

You can also find a contribution from Don Phillips, 
bringing actuarial math on the handheld together 
with some additions for the former DERIVE file. 

Peter Liicke-Rosendahl wrote a mail concerning 
some „gaps“ in my DNL# 15 contribution on Curve 
Discussions. We had an interesting exchange of 
mails and files and I'll come back to this issue in 
the next DNL. 

TIME GOES AFRICA. TIME 2008 will be held 
22 - 26 September 2008 at Buffelspoort Resort, 
South Africa. 

Michel Beaudin and I were asked to act as Co- 
Chairs for the DERIVE-Nspire-CAS-strand. Michel 
and I found it a good idea to have a completely new 
group of colleagues in the Program Committee. We 
would like to thank alle members of earlier PCs. 
They did a great job and they really deserve a 
break. 

You can find the first flyer and invitation (together 
with another announcement) at the end of this 
DNL. 

Finally I'd like to draw your attention to the pre- 
sented books (unfortunately only in German). If 
anybody of you knows about interesting publica- 
tions, then please let me know. 

By the way, Michel Beaudin provided a paper full 
of exercises for his students on the EMTS (in 
French). We will publish it in the next DNL. 

(Have a also look on SketchUp - page 20!) 

Best regards to you all 



Download all DNL-Derive- and Tl-files from 

http : //www . austromath . at/dug/ 
http : //www . bk- teachware . com/main . asp?session=37 505 9 
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The DERIVE-NEWSLETTER is the Bulletin 
of the DERIVE & CAS- 77 User Group. It is 
published at least four times a year with a 
contents of 40 pages minimum. The goals 
of the DNL are to enable the exchange of 
experiences made with DERIVE, 77-CAS 
and other CAS as well to create a group to 
discuss the possibilities of new methodical 
and didactical manners in teaching mathe- 
matics. 


Editor: Mag. Josef Bohm 
D'Lust 1, A-3042 Wurmla 
Austria 

Phone/FAX: 43-(0)2275/8207 
e-mail: nojo.boehm@pgv.at 


Contributions: 

Please send all contributions to the Editor. 
Non-English speakers are encouraged to 
write their contributions in English to rein- 
force the international touch of the DNL. It 
must be said, though, that non-English arti- 
cles will be warmly welcomed nonetheless. 
Your contributions will be edited but not 
assessed. By submitting articles the author 
gives his consent for reprinting it in the 
DNL. The more contributions you will 
send, the more lively and richer in contents 
the DERIVE & CAS-77 Newsletter will be. 

Next issue: December 2007 

Deadline 15 November 2007 


Preview: Contributions waiting to be published 

Some simulations of Random Experiments, J. Bohm, AUT, Lorenz Kopp, GER 

Wonderful World of Pedal Curves, J. Bohm 

Tools for 3D-Problems, P. Luke-Rosendahl, GER 

ANOVA with DERIVE & Tl, M. R. Phillips, USA 

Financial Mathematics 4, M. R. Phillips 

Hill-Encription, J. Bohm 

Farey Sequences on the Tl, M. Lesmes-Acosta, COL 
Simulating a Graphing Calculator in DERIVE, J. Bohm 
Henon & Co, J. Bohm 

Do you know this? Cabri & CAS on PC and Handheld, W. Wegscheider, AUT 

An Interesting Problem with a Triangle, Steiner Point, P. Luke-Rosendahl, GER 

Overcoming Branch & Bound by Simulation, J. Bohm, AUT 

Diophantine Polynomials, D. E. McDougall, Canada 

Graphics World, Currency Change, P. Charland, CAN 

Precise Recurring Decimal Notation, P. Schofield, UK 

Solving an Optimization Problem with the Tl-Nspire CAS, B. Kutzler, AUT 

Cubics, Quartics - interesting features, T. Koller & J. Bohm 

Logos of Companies as an Inspiration for Math Teaching 

Exciting Surfaces in the FAZ 

BooleanPlots.mth, P. Schofield, UK 

What is hiding in Dr. Pest? B. Grabinger, GER 

Truth Tables on the Tl, M. R. Phillips 

Advanced Regression Routines for the TIs, M. R. Phillips 

Directing Our Suspicions with AHP, C. & P. Leinbach, USA 

Where oh Where is IT? (GPS with CAS), C. & P. Leinbach, USA 

Groebner Bases, A. Perotti, IT 

Embroidery Patterns, H. Ludwig, GER 

Inspheres and Exspheres of a Tetrahedron, D. Sjostrand, SWE 
Complements de mathematiques (profil Genie Mecanique), M. Beaudin, CAN 


and others 


Impressum: 

Medieninhaber: DERIVE User Group, A-3042 Wurmla, D'Lust 1, AUSTRIA 
Richtung: Fachzeitschrift 
Herausgeber: Mag. Josef Bohm 









Conics Made Easy 

for TI-89/Titanium and TI-92+/Voyage 200 

DUG-Member Nils Hanhfeld offers two new packages for the Tl-handheld devices: CME 
(Conics Made Easy) and DME (Differential Equations Made Easy). Both products can be 
purchased via Internet: www.ti89.com . 

I'd like to present some screen shots from CME. I'll take the Voyage 200 because of the 
larger screen which needs a lot less scrolling the screens. 

After loading the Flash Application you are offered a bar with 5 drop down menus. The first 
menu gives the most important information about the conics, eg. the ellipse: 



The most powerful menu is F2 Analyze. You can enter either the coefficients of the equation 
of a given polynomial (in rectangular coordinates only) or the full equation (see the example). 



CME performs a complete analysis of the given conic followed by its plot. 
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CME switches automatically to DEG-Mode. 


Comparing the plot with the DERIVE plot (plotting the implicitly given expression) you will 
recognize that the Tl-plot has sometimes a scale adapted to the dimension of the conic and 
that the Tl-plot shows a nasty gap between the two parts, caused by problems plotting 
curves which are close to vertical runs. The ellipse is splitted into two parts. 



This problem could be overcome by plotting the full conic in parameter form. 

Here is a step by step approach for plotting the given conic in parameter form. (The last time 
when I performed a calculation like this was when I was 18 writing my end examination in 
gymnasium - several years ago ...) 


General fiorm ffx.y) = Ax a 2 + Bxy + Cy A 2 + Dx + Ey+F = 0 
2 2 

#1: f (x , y) := 4-x + 3«x«y + 5-y - 5«x + 3«y - 20 

First of all check the type! 

ch_comc(u, dis) := 

Prog 

dis := d(3(u, x), y) A 2 - 3(u, x, 2)«3(u, y, 2) 

If dis > 0 

#2: "hyperbola 11 

If dis < 0 
"ellipse 11 
"parabola" 

#3: ch_comc(f (x , y)) = ellipse 

Finding the center 


f 

f 

r d 

d 


S A 


" 59 

39 " 

cr := 

SOLUTIONS 

— f(x, 

y) = 0 a — f (x , y) = 0 , 

[X, y] 


= cr := 

i 



c 

^dx 

dy 

J . 



_ 71 

71 _ 






1J 





ft(x, y) := f(x + cr , y + cr) 

#5: 12 

This is the conic with its center in the shifted origin: 
(the equation does not contain any linear terms) 

2 2 
2S4-x + 213-x-y + 355-y - 1626 

#6: ftCx, y) := 

71 

We rotate the coordinate-axes that the axes of the conic 
will lie on the xy-axes , finding first the rotation angle ex 
from tan(2cx) = B/(A -C) : 
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ATAN 

#7: 

□ □ 

2 f(x. y) 

dy dx 


f 1 A 

ACOT 

L 3 J 

JV • — 

fd A2 fd yz 

— f(x, y) - - f(x, y) 

^ LdxJ Ld yj j 

LA ■ — 

L 2 J 

— IJ. 1 — 

2 


We apply the rotation -Formulae 

#3: f r(x, y) := ft(x-COS(oc) - y-5IN(ac)i x-SIN(oc) + y-C OS(oc)) 


#9: 


f r & ■ y) 


2 2 

71 «x -CjiO - 9) - 71 -y .(JiO + 9) + 3252 


142 



Finding the vertices by calculating the intersection points with the axes: 

#10: vertr(u) := APPEND(SOLUTIONS(u = 0 a y = 0, [x , y]), S0LUTI0NS(u = 0 a x = 0, [x , y])) 


#11: vertr(fr(x, y)) 


#12: vertr(fr(x, y)) 


2-7(313-710 + 7317) 
71 

2-y(S13.yiO + 7317) 
71 


0 


0 


0 

0 

2-^(7317 - S13-710) 
71 

2-7(7317 - 313-710) 
71 


2.301072673 

-2.301072673 

0 

0 


0 

0 

1.940611662 
-1.940611662 _ 


We can obtain the normal form using Eigenvalues, too. But this was not the way I did it in 
1963 in gymnasium. (I'll do this in the next DNL.) 
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The parameter form with coordinate axes shifted and rotated 
(& length of half axes): 

#13: [a := 2.S01072673, b := 1.940611662] 

#14: [a-COS(t), b-SIN(t)] 

To obtain the parameter form of the original conic we have to 
rotate and shift back: 


#13: 


#16: 


fp(t) := [a.COS(t), b-SIN(t)]. 


COS (-tx) - SIN(-a) 

SIN(-ot) COS(-ffl) 


fp(t) := 


2801072673-7(5 -JiO + 50)-COS(t) 970303831-^(30 - 3 -JiO) ■ SIN(t) 

+ + 


10000000000 
2801072673-^(30 - 3 -7l0) -COS(t) 
10000000000 


3000000000 

970305831-7(5-710 + SO)-SIN(t) 39 

5000000000 71 J 


59 


71 


Having plotted the conic in parameter form, we want to add the center and the vertices (and 
maybe for hyperbolas the asymptotes, too). You can also add the foci. 


#17: conpts := VECTOR 


u + cr, u, APPEND(vertr (f r(x , y)), [[0, 0]]) 


#18: conpts = 


#19: conpts = 


4(77235 -fLO + 406500) 59 

+ 

355 71 

59 7(^7235-710 + 406500) 


71 355 

7(406500 - 77235-710) 


59 


59 


71 


355 71 

7(406500 - 77235-710) 


355 


59 


71 


3.103334232 

-1.441362401 

1.965681513 

-0.3037096823 

0.8309859154 


-2.187111775 

1.088520225 

1.025010271 

-2.12360182 

-0.5492957746 


COS(-ffl) - SIN (-a) 
SIN(-ffl) COS (-a) 



When I worked the first time with CME I missed two options: 

- A general conic is given by five points: Enter the points eg in form of a matrix and receive 
the equation of the conic containing the points - followed by its analysis + plot. 

- All conics can be presented with one common function - in polar coordinates. This would 
be a welcome opinion to make students more accustomed with polar forms. 
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I suggested Nils to implement these features and he promised to do. I'd like to present how 
to find the equation of a conic which is given by five points. 

Starting with the general equation the problem is to solve a system of five linear equations 
with five unknowns - this should be no problem for a CAS: 

The screen shot shows 5 points (matrix m_), then I call function conic(m_) and receive im- 
mediately the equation of the conic - in rectangular form. Finally I plot its graph. 


rjSIz^Il^celRe^aphlM^hlD^r^l:^ 1 


□ 


□ 

□ 

MAIN RAD AUTO FUNC 


^IfluTbrafcLTc 

f FHt Y FE T FAt 1 

Other |PrgnI0|C lean Up 

LI 

■ zerosl. - . 0752565, 
{l.223-(jx 2 - 5. 1 
■{ 1.222592004584- 

■{l ,222592004584- 

36442426 "x^ + x ■(. . 1622 
07 -x + 6.994 - . 690 (x 

+ ► 

4 [Jx 2 - 5. 1 073280986667^ 
Done 

4 (jx 2 - 5. 1 073280986667^ 
Done 

ans<2M2]-Kp2<x>l 

MAIN RAD AUTD FUNC E/30 


1 gebrafca lclDtherlprgn I oic 1 ean UpT 1 

1 I 

2 -3 


1 I 

2 -3 

■ 0 3 

+ n_ 

0 3 

2.5 -1.5 


2.500 -1.500 

.6 "I . 


.6 "I 


■ conicO-} 


^5 ■ x 2 + x ■( . 162 ■ y + . 772) + . 096 ■ y 2 + . 04=> 

■■■5419ttv A 2+ .04484986697075*^-1 

MAIN F: H [i AUTD FUNC E/30 


["^Iz^Il^calRa^aphlM^thlD^I" ^ V \ 





MAIN RAD AUTD FUNC 


Function conic() is only a transfer of a DERIVE-function. See the function and an “over- 
painted” picture of a fountain using the derived equations: 

2 2 

#1: contx, y) := a-x + b-x-y + c-y + d-x + e-y = 1 


# 2 : 


com'c(pts) := (5QLUTIGNS(VECTQR(con(v , v ), v, 

1 2 


Pts) , [a, b, c, d, e])) 



2 


y 


j 


= 1 


1 


#3: com'c(elll) 



I was very happy with this function for a long while, but a “bug” is hidden. Do you see what is 
missing? No, then try to find the equation of a conic which is passing the origin! 


In the meanwhile we had a fruitful cooperation and many of the ideas have been 
realized for an updated version of CME. 
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In the next DNL I'll include a DERIVE-function (program) performing a complete analysis. 


2 2 " 
3 -5 


#33: pp := 0 0 


6 0 


L -3 -1 J 


#34: corric3?(pp) 


#35: 


hyperbola, x 


2 446 ■ x ■ y 


- 6-x 


3-y 1242 -y 


95 


19 


95 


S3 S3 


2 35 


3019 3019 


4 -73019.x 223 -x 


15 


15 


11134-73019 207 


4 -73019 -x 223 -x 11134-73019 207 


15095 


15 


15 


15095 


Information provided by Bernhard Kutzler: 

Hi all, 

Did you know that Derive has an entry in Wikipedia - with topical information about the dis- 
continuation and replacement with Tl-Nspire? 

Look at: http://en.wikipedia.org/wiki/Derive computer algebra system 
and in German: http://de.wikipedia.org/wiki/Derive 

Best regards, 
bernhard 


Restricted use of Slider Bars in DERIVE 

Heinrich Ludwig, Germany 

Preparing a visualisation in the frame of an introduction into Graph Theory for my students I 
wanted to turn on and off edges connecting nodes using slider bars. This works very easy 
presenting the edge as a segment connecting two nodes by a 2 by 2 -matrix and multiplying 
the matrix by a factor a with possible values 0 or 1 for a. With a = 0 the segment shrinks to 
length 0 and with a = 1 one can observe the full edge. 

Everything seems to work properly: one can introduce many slider bars, but ... only the first 
nine of them will be really working. I found no way to overcome this deficiency. 

DA/L: / wrote to Albert Rich and Theresa Shelby. The answer confirmed this restriction. It 
would need changing the DERIVE code plus recompiling. Unfortunately we cannot expect a 
new release, so we have to live with only 9 working slider bars. 


D 
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Ignacio Larrosa Canestro 

I know that Tl is not continuing the development to Derive, but I expect that this list will con- 
tinue alive. Someone may have some explanation for the slowness shown by Derive execut- 
ing this function: 

With argument k = 2 A m + 1 , for values of m between 1 and 15 delays it needs only fractions 
of a second. But from m = 15 it needs much more time: until 182 seconds for m = 19, and 
even to exhaust my patience for m = 20, on a Pentium IV 2800MHz. 

fflO := 

Pnog 
k := k/2 
m := 1 
Loop 

#1: If INTEGER?(k) 

RETURN m 

k := k-CEILING(k) 

DISPLAY(DIM (FLOOR (k))) 
m :+ 1 

Really alone it should do m iterations, the answer has to be m + 1 , as really gives for m less 
than or equal to 19. 

Thanks in advance, 

Ignacio Larrosa Canestro 
A Coruna (Espana) 

ilarrosa@mundo-r.com . , . 


Albert Rich 

Hello Ignacio, 

I don't think it strange at all that your function is so slow in computing f(2 A 20+1 ). Including 
the statement 

DISPLAY (DIM (FLOOR (k) ) ) , 

in the loop of your function displays the number of digits of k for each iteration of the loop. 

As you will see the number of digits doubles with each iteration. By iteration 20, multi-million 
digits numbers are being multiplied together! 

In fact, your function is actually a good example of the high performance of Derive's infinite 
precision arithmetic. 

Aloha from Hawaii, 

Albert D. Rich 
Co-author of Derive 

Peter Luke-Rosendahl, Germany 

Hello Josef, 

...I like to read your revised reprints. I am just trying your „Discussion of a Curve“ from revised 
DNL#15. I found out that “terrace points” (= points with f = f ’ = 0) are not recognized properly. The 
reason might by that within an IF-condition sign(l) and sign(O) are not distinguished . . . 

In his next mail Peter sent his proposal using the fact that the sign of the first derivative 
changes in a turning point, but does not change in an extremal value with slope = 0. 

I implemented Peter's idea and by the way I tried to improve the “program” from revised 
DNL15 to also consider points with vanishing higher derivatives (giving “flat points” ...) 

Many thanks, Peter for the fruitful discussion. 
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All bodies are falling equally fast - 
Are all bodies falling equally fast? 

Josef Lechner, Viehdorf, Austria 


Introduction 

More and more we are faced with the demand 
on real life applications in teaching mathemat- 
ics. But teaching real life applications has the 
consequence to deal with complex models. As 
one example we will investigate the free fall 
under considering air resistance [1], which is 
often discussed in school but only rarely really 
treated. 

Equation of motion 

In contrary to free fall in vacuum a body falling 
free in the air is not only influenced by the 
gravitational force but also by a retarding resis- 
tance force, which is described by Newton's 
formula for the resistance 


F w =c w A'- 


Equation of motion follows 
m • a = m • g - c w A- 


2m 

d t 2m 


The numerical approach 

Supported by a CAS (DERIVE, Voyage 200, 
TI-92+ or TI-Nspire) it is easy to obtain very 
quick a numerical solution. (In this paper we use 
the TI-92+ / Voyage 200, but all procedures can 
be reproduced using DERIVE or the Nspire.) 
When we enter the equation of motion in the 
Y=Editor (having set the DIFF EQUATIONS 
Mode) we receive the respective direction field 
in order to have a first overview. (Drag coeffi- 
cient c w = 0.8, cross-sectional area of the falling 
body A = 0.8 m, air density p = 1.3 kg • m" 3 and 
mass m = 70 kg: the resulting coefficient of the 
quadratic expression is 0.006.) 



Our initial condition is v(*=0) = 0. So we obtain 
the increase of velocity which - as expected - 
decreases and approaches a constant value. 



The end velocity which can be seen in the solu- 
tion curve will be reached when restistance of 
air equals gravity. Equating F g and F w yields the 
end velocity 



This is in our case 

v end ~ 40.63 m s~ x «14 6kmh~ l . 
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Rate of fall as a function of time 

For obtaining a presentation of the velocity 
function it is necessary to solve the respective 
differential equation. To have a better survey we 


set b := - 


2m 


-. For solving this DE the TI- 


devices provide the function desolve(), which 
makes possible to find exact solutions of most 
1 st and 2 nd order ODEs. 




j 

' 40.629067843678-3.6 


146.264644 
... 2 


v 2-n v v 2-n 

1 deSo 1 ue( v 1 = g-b-v^ and o£0} = 0 , i , o) | ► 
ln(|JE ■ o - Jg|) - ln(JE ■ v + Jg) 
2 ■ JE~g 


- = -t 


... and u<0> = 0 J . t J .u> I sf>0 and b>0 

MflIH FiflD flUTD DE E/30 


Using desolve() and defining g> 0 and b> 0 we 
receive the desired expression which can be 
solved for the velocity v. 


F£T T f::t 1 
| ▼ |fl 1 gebr a |Ca 1 c | 

r FHt Y FE Y FfiT Y ’l 

Other |PrgmI0|C lean Up| 

■ soloej 

a Jg- tai 

ln|JJb ■ u - J 

g| J - lnUb*v + JgJ _ x || 

L 2-JFi " "’ V J 

■ianh(Jb~iU Ji j. 

Jb =I,J f 2 jb^i t + 1 

nhQbglO . „ .8-.S-. 

|g-9.Sl and b- 140 r 

40. 629068 ■ ianhC . 241453 ■ t) 

13=9-8... 


DE IY3Q 


We can summarize the solution: 
v(0 = Jj- ■ tanh 


Back substitution under consideration of 



and 


'Jbg = 


g-c w A p _ g 


2m 


v end 


results finally in the velocity function 
v(t) = V end -tanh 


f \ 
V V end J 


( 3 ) 


Working with our data we can plot the exact 
solution (use the Draw = F6-tool) and observe 
an excellent match of both curves. 


_ Jg ■ tanhQb ■ g ■ l) , Jg. 

Jb 

. Jg ■ tanhQb ■ g ■ t) 


and 


2 ■ Jb-g ■ t + t 

. 8.8 


k 


-j= I 9 = 9-31 and b= 

40.629068 ■ tanh( . 241453 ■ t) 
■ DrawFunc 40. 629067843677 ■ tanh( . 241452> 

Done 

■ 3677*tanh< - 2 41 45 2 746 042 43 *x>l 

MAIN F;hD flUTD DE B/30 




^ 

> //? / // // // // // // // 
J ////////////////// 
/////////////////// 
/// / // // // // // // // // 
////////////////// 
I f/ ' / ' / ' / ' / ' ////// ' / ' / ' / ' / 777 ± 


FiliD flUTD 


DE 


Now let's have a look how DERIVE handles the problem: 
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( 1 . 3 - 0 . 8 - 0. 8 2 "I 

# 1 : DIRECTION_FIELD 9.81 y , x, 0 , 20 , 20 , y, 0 , 50 , 25 

L 140 J 


# 2 : SOLVE (DS 0 LVEl(g - b-y , - 1 , x, y, 0 , 0 ), y) 

2 -^/b-Jg-x 


2- Jb- Jg-x 

Jb- (e + 1 ) 



( ( 1 . 3 - 0 . 8 - 0. 8 

2 


# 4 : 

SOLVE DS 0 LVE 1 9.81 

L L i 4 o 

-■y i - 1 . ; 

/i o, 0 , y 

J J 


# 5 : y = 


0 . 48 29 05492 -x 

40 . 62906784 - (e - 1 ) 

0 . 482905492 -x 


( 0 . 482905492 -x "j 

40 . 62906784 - (e - 1 ) 

# 6 : TABLE , x, 0 , 22 , 0.05 

0 . 482905492 -x 

v. e +1 j 


It is interesting that DERIVE does not recognize the tanh in expression #3! We can read in 
the DERIVE Online Help: The hyperbolic functions simplify to equivalent expressions involv- 
ing exponentials. That is it! 


Way of the fall and its acceleration as func- calculator does not transform the expression into 


tions of time 

Now having derived the velocity function it is 
easy to find the distance of the fall (by integra- 
tion) and its acceleration (by differentiation). 

Integration of the function delivers the time- 
-distance function with a(£=0)=0. 

1^^— IfllgeSralc^^TotherTprgn I qTc 1 ean UpI 1 

■ |q (j| ■ tanhC-Jb ■ g ■ u))du | g > G and b>0 

l n (^ 2 -Jb“g t + l) - JtTg-t - ln(!2> 
b 

. InU 2 t -+ l) - JE~gi - ln(2) | L _ Q 

0 

ans<l> I t=Q 

MAIH DEG AUTD FAR 10/50 

Collecting the numerator in form of a log- 
function enables a better and more “beautiful” 
presentation of the desired function: 


s(t) = — • In 
b 


The calculator does not collect the logarithms by 
its own, but we can force it to do so: 


the cosh. We have to support the TI. 


)T Y F> T FH-f Y FF Y F St Y 

■ |fl 1 gebr a |Ca 1 c |Qther |Pr gn I Q |c 1 ean Up | 


-g- ■ ianhC-Jb ■ g ■ u) du | g > 0 and b 


< an s < 1 > I n= t > - < an s < 1 > I n= □ > 

MAIN DEG AUTD FAR 15/50 


)Y Y Y FH-f Y FF Y Ffi^ Y 

■ |fl 1 gebr a |Ca 1 c |Qther |Pr gn I Q |c 1 ean Up | 



■ 9 ■ U| 

f 2-Jb 

T^ u+ ! 



U4b 


ln [ f 2-.[b> 


b 


U 

2-Jb^-t +1 J| 



f 2jb T gt + 

2 e 4bTt V 

. 2 J 


i'T F^f YFJtY fh-f y ff t T 

■ |fl 1 gebr a |Ca 1 c |Qther |Pr gn I Q |c 1 ean Up 

2 .|^Jb-g- 1- 

"In — n . 


f 2 Jb - g ■ i 
b 


„ [ 2-^ b g t J ln(cQsh(JE~~g ■ i)) 


Back substitution using again v end gives 

2 f f M 

s(t) = ^^—- In cosh -^—-t 

v, , 


Let it find the indefinite integral, then substitute 
t and 0 for u and form the difference. The result 
is not the same as given in the formula and the 


g 


(4) 
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Graphic representation of our special example 
leads us to the typical run of a displacement 
function under consideration of air friction: the 
function is bent to left in the beginning and 
changes into a linear function (when the velocity 
becomes constant). 



Finally we get the acceleration function as the 
derivative of the velocity function. 



We perform back substitution using v end and 
obtain the formula for acceleration as a function 
of time 

a(t) = -P (5) 

cosh 2 — — • t 

V ^ end ) 

The displayed graph shows the expected run. 
Acceleration decreases very fast from the initial 
value g to zero. 



Of course, end acceleration a end and end velocity 
v end can be obtained analytically, too. 




9 .SI 


t-40.63 
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d 

# 30 : a(t) := — v(t) 

dt 

2 -g- t/v_ 

4 -g-e 

# 31 : a(t) :z 

2 -g- t/v_ 2 
(e + 1) 

2 ■ g ■ t/v_ 

4 -g-e 

# 32 : 

2 -g- t/v_ 2 
(e + 1) 

2 - 9 . 81 -t/ 40.63 

4 - 9 . 81 -e 

# 33 : 

2 - 9 . 81 -t/ 40.63 2 

(e + 1) 



#34: [g :e Real (0, «■) , b :e Real (0, «■) ] 

2- -/b- Jg- 1 

Jg-(e -1) Jg 

#35: lira = 

t-™ 2-^/b-^/g-t -/b 

'/b- (e + 1) 

li'm a(t) = 0 

t"KO 


# 36 : 
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Time of the fall as function of height of fall 

Determination of the inverse of the time- 
distance-function does not lead to the desired 
result. 



The reason for this incapability of solving the 
equation is the ambiguity of the arcosh(x): 



We have to support the CAS by deciding for the 
upper branch of the arcosh-function - negative 
times don't make any sense in connection with 
our model - then for x > 1 : 

arcosh(x) = ln|x + ^x 2 -lj. 

So we get from 


cosh (yfbg • t) = e bs 

Jbg ■t = \n[e bs 



finally 



And after introducing v end we obtain the formula 
for the fall(=drop) time 



We will compare the graphs of the functions for 
fall time with and without air resistance for 0 to 
100 m. 



Reference 

[1] G. Hepperger, Der freie Fall unter Einwir- 
kung des Luftwiderstands, Wissenschaftli- 
che Nachrichten, September 1991, BMUK, 
Wien. 


Josef received a letter from a German teacher, 
Madeleine: 

Completion: “Madeleine 's Problem ” 

Two parachutists, Klaus and Egon , leave the 
plane in an altitude of 2000m and they open 
their parachute in an altitude of 200m above the 
ground. Their drag coefficient c w is 0.45, the 
density of air (1.29 kg/m 3 ) decreases at a rate of 
5% per 100 m difference of altitude. Klaus 
keeps his body vertical and his acting cross sec- 
tional area of his body is 950 cm 2 . Egon is lying 
horizontally on the air bed, spreads his arms and 
legs reaching a cross sectional area of 7500 cm 2 . 
They both have a weight of 80 kg. Calculate the 
time which they need to reach the the moment of 
opening their parachutes and give reasons for 
your procedure. 


There are two forces acting on both parachutists: 

p • v 2 

F c =m-g (Weight) and F w =c w -A 

We have to consider that air density increases 
approaching the ground (g and c w depend also 
on the altitude, but this shall be neglected here). 

Leaving the plane in an altitude of H= 2000 m 
we can describe the density p as a function of 
the fall distance s: 

2000 -s 

p(s) = 1.29 -0.95 100 . 
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The resulting equation of motion is 


m ■ a = m ■ g - c w ■ A 


p 0 - 0.95 100 -v 2 


(H is the altitude of the jump). This leads to the 
following differential equation: 

c . A H MO 

^(0 + -r Po'0.95 100 -s'(t) 2 -g = 0. 

2m 

For our - Madeleine's - problem two differen- 
tial equations have to be solved - one for Klaus 
and another one for Egon : 


Klaus: 


0 . 45 - 0.095 


-• 1 . 29 - 0.95 100 —s\t) - 9.81 = 0 


0 45-0 75 2000-sco 

s\t) + : 1 . 29 • 0.95 100 - s\tf - 9.81 = 0 

2-80 

The initial conditions are s(0)=0 and v(0)=0. 


The DEs can be solved numerically (eg by ap- 
plying the Runge-Kutta-method. For editing 
these DE's on the TIs one has to take into ac- 
count some conventions: 

yl corresponds with s K [ aus (1 st DE) 
y2 corresponds with v K i aus (1 st DE) 
y3 corresponds with s Egon (2 nd DE) 
y4 corresponds with s Egon (2 nd DE) 

After entering 5 ’ and v’ 



we can read off the requested times on the 
graphic screen - at least approximatively. 




We can read off a time of 20.2 seconds for 
Klaus and 28.1 seconds for Egon, yl and y3 are 
the values of their falling distance. It is obvious 
that Egon needs more time 

So the problem is solved - at least numerically. 

It is interesting to investigate the time - velocity 
diagrams of both jumpers: 



We can recognize the braking influence of the 
increase of the air density during fall time. 

Finally we will have a look at the distance- 
velocity diagrams. At first it seems that we can- 
not have both diagrams in one graph window: 



Acitvate yl ’ and y3’ in the Y=-Editor. 

We can get a more accurate estimation using the 
table: 


MAIN 


RAD liUTD 









D 
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All these graphs are easy to create by using the 
Axes-Menu (F7) in the Y=-Editor. In Axes 
CUSTOM you can assign variables to the axes 
as you like. (The screen shot demonstrates how 
to set F7 Axes > CUSTOM in order to obtain 
the distance-velocity diagram for Klaus. 


Plot ■ X=5irv^UlJti\5tr1 y:5i|vni1flti\XP-rOZl 


Plot 

f HUES 


t0=c 

yl'=L 

yil=G 

fixes . . CUSTOM-* 

X fixis yl-* 

V fixis GUI* 

[ CEnter-=SfiUE") CESC=CfiNCELl>J 


4 



] 


yi2=0 


v2‘ q>=-.45tt.095/160ttl.29tt<.9... 

USE * AND ± TD OPEN CHOICES 


The differential equations and DERIVE: 


And with a little trick we can display both 
graphs on the same axes: Save both curves as 
separate pictures (FI -Menu > Save Copy as) 
and load both into an empty graph window. 



We follow the advice given in the Online Help how to treat second order differential equations 
in order to obtain a numerical solution. 

Expressions #1 7 and #1 8 deliver a table similar to the Tl-table (increment = 0.1). We can 
extract the coordinates for the desired solution curves or at first extract some rows of the 
table to have an estimation for the time needed until opening the parachute (after 1800 m): 


#17: klaus RK 


v, 9. SI 


0.45-0.095 


160 


(2000 - s)/100 2' 

-1.29-0.95 -v 


[t, 


v], [0, 0, 0], 0.1, 



#1S : agon 


RK 


v, 9. SI 


0.45-0.75 (2000 - s)/100 2 

1.29-0.95 -v 


160 


[t, 


v], [0, 0, 0], 0.1, 



The first column shows the time t, the sec- 
ond one the distance s and the third one 
gives the velocity v. 

We plot time the time-distance diagram, then 
the time-velocity curve and finally the veloc- 
ity-distance graph by selecting the respec- 
tive columns of the matrices klaus and egon. 


#23: 

klaus COL 

[i, 

2] 

#24: 

agon COL 

[i, 

2] 

#25: 

IS 00 




SELECT (17 50 < v < 1850, v, klaus) 
#19 : 2 


' 20 

1763 .940123 

153. 36S 7157 ' 

20.1 

1779 . 2S9996 

153.6273766 

20.2 

1794.665328 

153.8779058 

20.3 

1S10. 065305 

154.1202846 

20.4 

1S25.4S9112 

154.3544954 

. 20.5 

1S40. 935932 

154.5805225 . 


SELECT (1790 < v < 1810, v, egon) 
# 21 : 2 



pl8 


Josef Lechner: All bodies are falling equally fast? 


D-N-L#67 


#26: klaus COL [1, 3] 

#27: agon COL [1, 3] 


#28 : klaus COL [2, 3] 

#29: agon COL [2, 3] 




Endexamination 2001, Gymnasium Amstetten, Lower Austria 


a) 


b) 


The acceleration acting on a parachutist is equal to the acceleration due to gravity minus the re- 
verse acting brake acceleration caused by the air, which is proportional to the square of the fall 
velocity (prop, factor = 0.006). 

Set up the respective differential equation and give its solution as v(t) with v(0) = 0. 

(g « lOm/s; Work in Exact Mode!) 


Find the respective time-distance function with s(0) = 0: How far did the parachutist fall until the 
parachute opened, if he stayed 5 seconds in free fall? Which is his fall velocity at the moment 
when the parachute opens? 


m 1 ln(cosh (tJb-g)) 

Show that s(t) = - 

b 


also describes the distance falling down. 


c) Which end velocity could be reached? Using the ^-function find an expression for the end veloc- 
ity. 

d) Describe the run of velocity applying a restricted growth function. Choose an appropriate prop, 
factor that the function matches the velocity function from above as accurate as possible. 

e) Extra Credit: Increasing altitude is accompanied by a decreasing air density, which even influ- 
ences parachuting. Leaving the airplane in an altitude of 2000m above ground this influence can 
be modelled by 

2000 - 5 ' 

b(s) = 0.003 -0.95 100 . 

Solve the resulting DE approximatetively in the DE-Mode and give a sketch of the resulting ve- 
locity function. 


Remarks: This example uses the many features of the TI-92/V 200 in a high degree. The idea for this example 
gave a question posed by a certain Madeleine (German teacher) some time ago. She asked if a parachutist might 
become slower during his/her dropping down to earth. As part e) shows, he/she might because of the braking 
effect of increasing air density. 

Part a) cocnsits of setting and (exact) solving a DE. This could be done without CAS, but then we had to leave 
other parts of the example. 

Parts b) and c) aks at one hand for several basic Calculus techniques, and at the other hand students have to 
demonstrate their competence inshowing the equivalence of not elementary expressions. This compentence is 
very important in the frame of a CAS supported math education. 

In part d) the students shall show their capability of modelling and finally 
part e) makes use of the TI's capability of solving DEs approximatively. 
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A Repository of Mathematical Knowledge 

Albert D. Rich, Applied Logician 
17 August 2007 

Having not heard from me in a while, Josef recently asked what I have been up to since 
Texas Instruments terminated development of Derive. The following is my response: 

No longer having to worry about customer expectations and upward compatibility has had a 
tremendously liberating influence on me. I am now free to follow the math wherever it leads, 
and it leads to some pretty interesting and unexpected places. 

For example, it is now clear to me that a new arbitrary-element-of operator, along with open 
and closed intervals, is essential for representing the infinite sets that can arise when solving 
equations or simplifying Boolean expressions. This in turn leads to the need to algebraically 
simplify expressions involving infinite sets (including intervals) and the arbitrary-element-of 
operator. Since set theory provides the foundation for all mathematics, computer algebra 
systems (CAS) need to spend as much effort simplifying sets as they have traditionally done 
with scalars and vectors. 

I am in the process of collecting and categorizing the elemental transformation rules (mathe- 
matical factoids) necessary to intelligently automate the simplification of mathematical ex- 
pressions. Each rule must not only include the domain over which it is valid, but the circum- 
stances under which it actually simplifies an expression. Otherwise, application of the rule 
can make things worse or even lead to an infinite regress. This issue is more fully described 
in my talk “Automating the Simplification of Mathematical Expressions” delivered at the Visit- 
ME 2002 conference held in Vienna, Austria [1] . 

The knowledge in current computer algebra systems is hard-coded in conventional pro- 
gramming languages. To add new knowledge to such systems, a mathematician first has to 
teach a programmer the required math. Then the programmer has to translate the math into 
computer code. Each step in this process is time-consuming and fraught with the possibility 
of error. 

Instead, I propose building a repository of mathematical knowledge, the contents of which 
will be open and accessible to all. Rather than code, the knowledge will be abstracted in the 
form of transformation rules stored in a tree-structured discrimination net. Mathematicians 
and scientists will be able to add new rules to the net, generalize existing ones, and fix or 
delete invalid ones. 

This is analogous to the way users add knowledge to Wikipedia, the popular on-line encyclo- 
pedia. But to solve problems using Wikipedia, or any reference manual for that matter, the 
reader must manually find and then apply the required knowledge. In the repository I envi- 
sion, the knowledge is in a vastly more useful form since the rules appropriate for particular 
problems are automatically found and applied. 
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I have no desire to create and sell yet another system in the already crowded CAS market. 
Rather, my intent is to prove the viability of using a discrimination net to store mathematical 
knowledge. To provide a proof-of-concept of this paradigm shift, the following components 
needs to be built: 

1 . A "starter kit" of transformation rules sufficient for doing mathematics up to college level 
calculus and vector algebra. This will probably require 5000 to 6000 rules. Thus far I 
have collected about 4300. The rules are expressed in standard math notation as linear 
strings of text. 

2. A computer program that can transform these text strings into parse trees, and then build 
a discrimination net by inserting the parse trees at the appropriate place in the net. 

3. A computer program that can simplify mathematical expressions by repeatedly finding 
and applying applicable rules in the discrimination net until no more rules apply. 

Items 2 and 3 will be a challenging, but relatively straight-forward, programming project since 
all the mathematical knowledge is encapsulated in item 1 . They can be written in any mod- 
ern programming language (e.g. Java) that provides automatic garbage collection of linked- 
list data structures and libraries for infinite precision arithmetic. Using such a language has 
the advantage of making the repository platform independent. Since the rules are stored in a 
discrimination net, the appropriate rule to apply can be quickly found in log n time, where n is 
the number of rules in the system. So efficiency should not be a problem so long as the 
rules are carefully crafted. 

Once the advantages of such a repository become recognized, I am convinced it will be 
adopted as the 21st century standard for storing mathematical and scientific knowledge in a 
readily accessible form. Hopefully one or more DUG members will help me make this ambi- 
tious project a reality. 


[11 You can find this lecture among the files which accompany this DNL. Thanks Bernhard Kutzler for 
permission to include this paper. 

As I wrote in my Letter of the Editor, it would be great if DUG members would help realising Albert's 
ambituous plans. Josef 


Do you know Google's SketchUp, which is a 
great program for modelling 3D objects. 
SketchUp can be downloaded for free from 
Google's websites. One of my students, Dr. 
Erwin Rybin, delivered an excellent 
SketchUp tutorial and permitted its publica- 
tion. You are invited to download this tuto- 
rial. It is a pdf file written in German, but it 
contains so many pictures that you might 
understand it even with little or no knowl- 
edge of German. 



Pantheon in Rome, rebuilt with SketchUp. 
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Don Phillips provided an extended contribution on Actuarial Math with DERIVE 6 in 
DNL#57. Here is the version for the handheld. Don sent a grouped file for the 
TI-89/TI-Titanium. It was no problem to install the files on the TI-89/Titanium and on the 
TI-92+/Voyage200 as well. In Don's paper there are only TI-89 screens - according to 
the title of the contribution. The TI-92-screens might be a bit better readable. I don’t want 
to print the whole paper, because there are many similarities with the DNL#57 paper but 
it will be among the zipped files accompanying this DNL. I'd like to give a selection of ex- 
amples comparing the DERIVE procedures with the Tl-procedures. 

When Don sent the paper he wrote that he “added a few extra things that weren't in his 
Derive file”. I tried to reproduce these “extra things” for DERIVE. Let's start with Don's in- 
troduction. 


Actuarial Math on the TI-89/TI-92+/Voyage 200 - 
- and DERIVE as well 

MacDonald R. Phillips 
don.phillips@gmail.com 
July, 2006 

The TI group file, actuary.tig, provides a set of basic functions, such as D, N, S, C, M, and R, for 
computing the actuarial present values of life insurance and life annuities, as well as premiums and 
contributions. It assumes a basic working knowledge of actuarial mathematics. The calculations here 
represent net, not gross, premiums and contributions; i.e. they do not include the expenses of the in- 
surance company or agent commissions, etc. (The complete TI-89 paper is among the files). 


To set up the custom menu for the actuarial functions, 
change to the ACTUARY folder and run the MENU pro- 
gram. Press FI, 1 to access the CommFunc(). 


Enter the assumptions: 


^ 0s 

ANNUIT I EsJaDU E ANnI^TI 



2!Extrapol O 
3: CurInfo< ) 
4: qxmy 
5: y. 

6 : ID< 

7: ClrHone 
8: Custnoff 



■ nenuQ Done 


menuO 


TYRE DR USE * [ENTER]=DK AND [ESC]=CANCEL 


As noted, the mortality tables and age setbacks are entered 
in a 2-row matrix. Be sure to enclose the mortality table 
names in “ The names of the four mortality tables in- 
cluded with the program are located under F2. You can 
use them to enter the names in the matrix. If you add other 
mortality tables to your calculator, you can edit the MENU 
program to include their names under F2. 

The a2000m and a2000f tables are the annuity tables #864 
and #863 in the Society of Actuaries Table Manager Soft- 
ware. 


f FiT T F£t Y Y fit Y FEt 

|Tools|Tabl 0 s|Ps & Qs|aNNUITIES|aDU ANN!* | | 

■ menuO 

Done 

IT M nupl994" 
■comnfunc 1. 

LL f up 1994 

6, 2, to, 10 9 J 

Done 

... .0; "f up1994" .0] .6 .2 . toj_10 A 9> 

ACTUARY RAD AUTD 

FUNC 1/40 


IPs & Qs ANNUITIES ADU ANN i 


2: a20 

3:mupl994 

4:fupl994 


It takes a good three to five minutes to calculate all the 
commutation functions for X, Y, and J. As the program is 
running, it displays ‘Processing x. . . etc., until it is done. 


■ rcenuQ Done 


menuO 


TYRE DR USE * E E H T E F; ] = D K AND [ESC]=CANCEL 
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The other functions under FI: Tools are Extrapol(), CurInfo(), and ID(i, nth, mth). The Extrapol() 
lets you choose between a uniform distribution (udd) and a linear distribution (lin) of deaths between 
integral ages. 


CurInfo() displays a list of the current assumptions. Press ENTER when you are done reviewing the 
assumptions. 


ID(i, nth, mth) converts between interest rates or discount rates compounded nthly to rates com- 
pounded mthly. For instance, a 6% annual effective rate of interest is equivalent to 5.84106. . .% com- 
pounded monthly. 


Finally, qxmy returns the age difference between x and y and % is simply the percent function. I just 
like having it in a convenient place and Custmoff restores the menu bar. 


[»lvK 

-1 

: 1 ■ ' ,[ 1 

r 11 roup 1994 11 

G 1 


L 11 f up 1994 11 

oj 


Int. Rate 

= 6 


flge Dif = 

2 


To:' = to:' 



Radix = 1GGGGGGGGG 


ACTUARY 

RAD AUTD 

FUNC E.-'Hii IXTH? 


( fi-tT f£t t f:=t t fh-t r fs- ire T ) 

| T 1 = I T .3 1. 1 s; S I F- S & T! = 1 H N N U 1 T 1 E L" I H [1 hNN|P | 1 

■ id(6, 1 , 12) 

5. 8410606784 

■ id<6, 1,«) 

5.3268903124 

■ id(6, 12, -12) 

5.9701492537 

HiMhPwJ 

ACTUARY DEG AUTD 

FUNC 3/30 


A 6% annual effective rate of interest is equivalent to 5.84...% compounded monthly, to 5.82...% 
compounded continuously and to a 5.97. . .% discount rate compounded monthly. 


1. Life Annuities 

There are 10 commutation functions for basic annuities. 
They all take at least two arguments: age and person. Per- 
son is entered as x, y, or j and is always the last argument 
entered. The three basic annuity commutation functions 
are Dx, Nx, and S. These are used to compute annual an- 
nuities. 


Dxm, Nxm, and Sm compute annuities payable mthly, e.g., monthly, quarterly, etc. sDx and sNx are 
used to compute annuities payable yearly that increase by a given percentage each year; the s stands 
for a salary scale. And sDxm and sNxm are used to compute annuities payable mthly that increase 
each year by a given percent. 


[t ool s]t ab lesjps Tqs] 

jsETHSraHi 

2! Dxm 
3:sDxC 
4: sDxnC 
5:Nx< 

6: NxmC 
7: sNx£ 

8: sNxm< 
9:SC 

IflDU^flNNftT] 




TYPE DR USE 4+ti + [ENTER]=DK AND [ESC]=CANCEL 


Example 5: Suppose X, age 25, wants to provide for a retirement annuity of $2,500 a month when he 
retires at age 65. What is the APV {Actuarial Present Value ) of this deferred annuity? If X wanted to 


2500-12- NM (65 , 12) 

#34: = 25784 , 22814 

D(25) 

2500-12- NM (65 , 12) 

#35 : = 1641.925657 

N(25) - N(65) 

2B0O-12- NM (65 , 12) 

#36: = 1687.477924 

NM(25 , 12) - NM (65 , 12) 


1 gebralca^lclotherlprgm I ole 1 Up! 1 

250G- 12nxm(65, 12, x) 

25784.23 

1641.93 

1687.48 

dx(25,x) 

B 25GG- 12nxn(65, 12, x) 

nx(25 , x) - nx( 65 , x) 
n 25GG- 12-nxm(65, 12, x) 

nxro(25 , 12 , x) - nxro(65 , 12 , x) 

.„/<nxm<25,12 J .x}-nxin<65, 

12 -x>>l 

ACTUARY RAD AUTD FUNC 3/40 


Thus, to receive $2,500 a month for life when he reaches age 65, X would only have to pay 
$25,784.23 at age 25. X’s contribution to his retirement fund would be $1,641.93 payable at the be- 
ginning of each year until he was age 64. Or, X could make an annual contribution of $1,687.48 with 
1/1 2 th ($140.62) payable at the beginning of each month until he was age 64 and 1 1 months. 
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Now, suppose X wanted to make annual contributions that 
increased by 10% a year to his deferred annuity. Here we 
have to use the commutation functions that incorporate a 
salary scale. 

X’s contributions begin at $298.84 when he’s 25 and end 
at $13,525.33 when he’s 64. 

Let’s take this calculation apart. First, the 2500*12 is the annual pension or annuity. Second, 
Nxm(65,12,x)/Dx(25,x) is the APV at age 25 of the pension. Next, (sNx(25,10,x) - 
- sNx(65,10,x))/sDx(25,10,x) is the APV at age 25 of a 40-year term annuity of 1 that increases by 
10% a year. Divide the APV of the pension by the APV of the annuity of 1 gives the beginning an- 
nual contribution to fund the pension which will increase by 10% each year following the first year. 


|Ti7*riT FT? Y f::t y fht y F? t f?? Y "| 
I t 4 -— |Fl 1 gebra |Ca 1 c |other |Prgm 1 0 |c 1 ean Up | | 

1 1 i- - 1 ? 1 1 


2500-12-nxn(65, 12, x) 
nxn(25 , 12 , x) - nxn(65 , 12 , x} 
2500- 12-nxn(65, 12, x) 

dx(25, x) 

snx{25 , 10 , x^ - snx(65 , IQ , x) 
sdx(25 ? 10, x) 

40 


293. S41 12019005 ■( 1,1) 

fB5t 

I ACTUARY 


ans<l>*l.l A 40 


1637.43 


293.34 


13525.33 


This last problem was not posed in DNL#57 and both functions snx() and sdx() don't have 
their equivalents in the DERIVE utility file. I tried to add these two functions making it easier 
for me by assuming nontruncated tables. 


#40: SD(age, sc, tab := x) 
#41: SN(age, sc, tab := x) 


^ sc ^|age 

1 + ■ D(age, 

L 100 J 

DIM (tab) - 1 

1 SD(i , sc 


tab) 

tab) 


i = ag e 


2300-12- NM(6B , 12) 

D(2S) 

#42: = 29S.S4H201 

sn(2B , 10) - sn(6B , 10) 

sd(2B , 10) 

40 

#43: 298. 84112 GO. -1.1 = 13B2B. 32662 


2. Advanced Annuities (all following examples are from the Handheld-Paper) 


The advanced annuity functions are located under 

F5: ADV ANN. 


Example 7: nEx(age, n, person) is the endowment func- 
tion; it is the APV of 1 to be received in n years. What is 
the APV of a $10,000 endowment to be received by X age 
35 in 5 years, if he’s still alive? 

nEx is not implemented in the DERIVE version, we 
will do this! 


( F1t Y Y FI:t Y FHt ^ 

|Too 1 s |T ab 1 es |Ps & Qs|flHNUITIES 

EaBEaufr 

0 


2: CaM{ 
3:CLM< 
4:CL12< 
5: JSMC 
6: JS12t 
7: sLM< 
3:sL12C 



TYPE DR USE + [ENTER]=DK AND [ESC]=CANCEL 


|ti7»T|T FEt Y FjT Y FHt Y FS Y FbT T "| 

| t f-— lo 1 gebra |Ca 1 c |other |Prgm 1 0 |c 1 ean Up | 

■ 10000 ■ nex(35 ? 5 ? x) 

7436. 15 

10000#nex<35,5,x> 

ACTUARY RAD AUTD 

FUNC 1/H0 
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New Function NEX 
Example 7: 

D(age + n, tab) 

#74: NEXfage, n, tab := x) := 

D(age, tab) 

#75: 10000- NEX(35, 5) = 7436.149361 

Example 8: 

#76: 1500 -CaM(5 , 12) = 6522.070426 

#77: 1500 -CaM(5 , -12) = 6490.477699 


Example 8: CaM(n, mth) is the present value of an annuity due certain of 1 payable mthly. If mth is 
negative, it is an annuity immediate, not an annuity due. Since this is an annuity certain, it is good for 
all ages and persons. What is the present value of an annuity due of $1,500 payable monthly for 5 
years? The annual interest rate, of course, is 6%. 


If this was an annuity immediate, the present value 
would be 6490.48 

These values can be computed with the Finance App, 
too: 



■ 15GG-caK5, 12) 6522. G7 

■ 15GG ■ caK5 ? - 12) 649G. 48 



ACTUARY RAD AUTD FUNC E/HO 


The annual annuity of $1,500 is payable monthly which is $ 125/month. It is received for 5 years or 60 
months. The payments are made at the beginning of the month and the 6% interest rate is com- 
pounded annually (left) and if the payments are made at the and of the month (right): 


It oo 1 sjconputel 1 


( F1t Y F£ Y ) 

|Tools|Conpute| | 

H=60.000 


N=60.000 

1^=6.000 


1^=6.000 

BFU=KdMKHBl 


BFU=B£KMCf^ 

PMT=125. 000 


PMT=125. 000 

FU=0.000 


FU=0.000 

PpY=12.O00 


PpY=12. 000 

CpY=1.000 


CpY=1.000 

PMT: END IrirffllBI 


PMT:@IS BEGIN 

Ft hint ualut 


uolut 


Instead of Finance App we can use TVM for Derive 6 from DNL#55: 


7VM(60, 6, pv, 125, 0, 12, 1, 1) = 


N 

IR(S0 

*PV 

PMT 

FV 

PpY 

CpY 

Payment 


60 

6 

-6522 . 07042711008 
125 
0 

12 

1 

BEGIN 


D 
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TVHC 60 , 6 , pv, 125 , 0 , 12 , 1 , 0 ) 


N 60 

IR( 9 Q 6 

*PV - 6490.47769871834 

PMT 125 

FV 0 


PpY 12 

CpY 1 


Payment END 


Example 9: CLM(age, n, mth, person) computes the 
APV of an n-year certain-and-life annuity of 1 pay- 
able mthly for a person age x. What is the APV of a 
10-year certain-and-life annuity of $25,000 payable 
quarterly for Y age 65? (In this case Y receives the 
$25,000 for 10 years, alive or dead, and for the rest 
of her life after 10 years if alive.) 

Example 10: CL12(age, n, person) computes the 

APV of an n-year certain-and-life annuity of 1 pay- 
able monthly for a person age x. If the above annu- 
ity were payable monthly instead of quarterly, the 
APV would be $292,729.67. 



■ 1500-caKS, 12) 6522.07 

■ 15G0 -caK5, -12) 6490. 48 

■ 25000 -clK65, 10, 4, y) 294659.31 


|25000ttclnK65,10,4,y> 


ACTUARY RAD AUTD FUNC 3/HQ 


The APV is $294 659.81. 


? i^Tn 1 g & b r ale a f =jo the r|p r g n I olci ear. Ur 


u 


' 1500 -caK5, 12) 

' 15Q0 -caK5, -12) 
'25000-C1K65, 10, 4, y) 
' 25000- cl 12(65, 10, y) 


6522.07 

6490.43 

294659.31 

292729.67 


25000ttcll2<65 J 10 J y> 


Example 11: JSM(age, percent, mth) computes the APV of a joint-and-survivor annuity payable 
mthly. It assumes that the annuity belongs to person X and person Y receives a survivor annuity if X 
dies. Age is the age of X, percent is the percentage of X’s annuity Y receives if X dies, and mth is 
how often the annuity is paid each year. For example, X age 60 retires and receives a pension of 
$65,000 a year payable monthly at the end of each month. X’s partner Y receives 55% of X’s pension 
if X dies. X is 2 years older than Y. What is the APV of the pension? 

The APV is $751,319.59. Notice that mth is nega- 
tive. This is because the annuity is paid at the end of 
the month instead of the beginning of the month. 

JS12(age, percent) computes the APV of a joint-and- 
survivor annuity payable mthly at the beginning of 
the month. What is the APV of the above pension? 

The APV is $756,736.26, over 5,000 more than if it was paid at the end of the month. This is reason- 
able since X gets his pension sooner rather than later. 


itiTUTlT F£T Y FjT Y FHt Y FE Y Fb^ Y ) 

f— - IfllgebralCalclotherjPrgmlolciean Up| 1 

■ 6500G- jsn(65, 55, -12) 


751319.59 

■ 65000- jsl2(65, 55) 


756736.26 

65OO0*.isl2<65.55>l 

ACTUARY RAD AUTD 

FUNC E/HO 


25000 -CLM( 65 , 10 , 4 , y) = 294659 . & 1 B 0 
25000 - CL 12 ( 65 , 10 , y) = 292729 . 66 S 5 
65000 - ]SM ( 65 , 55 , - 12 ) = 751315.5914 
65000 - ]S 12 ( 65 , 55 ) = 756736 . 2 B 7 & 
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Example 12: sLM(age, mth, scale, person) computes the APV of a life annuity that increases each 
year by a scale factor and is payable mthly. Y is age 55 and purchases a life annuity of $1,000 that 
increases by 2.5 percent a year and is payable at the end of each month. What is the APV of this annu- 
ity? 

The APV is $17,758.35. If there was no scale factor, the APV would be $13,406.23. This can be 
confirmed with the basic commutation functions. 

If this annuity was paid monthly at the beginning of each month, the function sL12(age, scale, person) 
could be used to compute it. 

As usual, the APV is a little higher because the payments come at the beginning instead of the end of 
the month. 

SLM and SL12 are new! 

1000- 5LM(55 , -12, 2.5, y) = 1775S. 35241 

1000- 5LM(55 , -12, 0, y) = 13406. 22575 

1000- NM (55, -12, y) 

= 13406.22575 

D(5B, y) 

1000- 5L12(5B , 2.5, y) = 17S77. 33715 


fFi-TWjY Y FjT Y FHt Y FE Y Fb^ Y 1 

| t |fl 1 gebra |Ca 1 c |other |Prgm 1 0 |c 1 ean Up | 

■ 1000 ■ sln(55 , -12, 2.5, y) 


17758.35 

■ 10O0-slm(55, -12, 0, y) 


13406.23 

1000 ■ nxn(55 , "12, y) 
dx(55 , y) 


13406.23 

■ 1000 si 12(55, 2.5, y} 


17877.34 

1000#sll2<55 -2 -5 

ACTUARY RAD AUTD FUNC 

H/HO 



3. Pensions 

Under F6 are some basic functions for computing 
the normal cost and actuarial liability of basic pen- 
sion. 

(These functions are not included in the actuar- 
ial file from DNL#57. Some of them make also 
use of the TVM-file from DNL#55.) 


The normal cost (NC) of a pension for any given year is the actuarial value of the part of the total pen- 
sion benefit assigned to the year following the valuation date, assuming valuation at the beginning of 
the year. The actuarial liability (AL) is the current value of past normal costs. 

Example 13: Traditional Unit Credit (TUC) 

The TUC cost method is usually used with pension plans that provide a flat pension benefit, such as 
$50 a month per year of service. The normal cost of this type plan is computed with the NCUC(bx, 
mth, cage, rage, person), ‘bx’ is the annual pension benefit earned in the following year, ‘mth’ is how 
often the pension is paid, ‘cage’ is the persons current age, ‘rage’ is the expected retirement age, and 
‘person’ is X, Y, or J. 

X earns a pension benefit of $75 a month for each year of service. If X age 35 expects to retire at age 
65, what is the NC of his pension in the following year? What is the pension plan’s actuarial liability 
if X was hired at age 30? 


t ihsT 


■ menuQ Done 

menuO 

TYPE DR USE + EEHTER]=DK liHD [ESC]=CAMCEL 



3: RLUC( 
4:NCEflNLD< 
5:fll_EFlNLD{ 
6: HCEflHLPt 
7:flLEFlNLF< 


D 
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The NC is $1,396.91. Notice that this is the amount 
that must be put into the pension plan at the begin- 
ning of the year to fund the benefit earned during the 
year. Since the benefit is yet to be earned, there is no 
liability associated with it until the end of the year. 


To find the actuarial liability use the 

ALUC(bx, mth, eage, cage, rage, person) function 
where 6 eage’ is the entry age into the pension plan. 


The AL is $6,984.54. This is 5 times the NC for X’s 
6 th year of employment. There is no liability for the 
6 th year because it has yet to be earned. 


And in DERIVE: 


|?Tpehs ionsTl i fe'i nsT 

1 

■ ncuc£75 ■ 12 , 12, 35, 65, >0 

1396.91 

NCU(X75#12 .12. 35. 65. x> 


ACTUARY RAD AUTD FUNC i/40 



|?Tpehs ionsTl i fe'i nsT 

1 

■ ncuc(75- 12, 12, 35, 65, -a) 

1396.91 

■ aluc(75 ■ 12 , 12 , 3G , 35 , 65 , >0 

6934.54 

ALU(X75#12,12,30,35,65,x 

B 

1 ACTUARY RAD AUTD FUNC 1/40 I 


|?TpENS IONsTl I FE ' T NS 

1 

■ ncuc(75- 12, 12, 35, 65, >0 

1396.91 

■ aluc(75 ■ 12 , 12 , 3G , 35 , 65 , >0 

6934.54 

■ 1396. 9G72677S93 ■ 5 

6934.54 

11 3 9 6 ■ 9 072 6 778 9 3 *5 Ml Mill H III 


1 ACTUARY RAD AUTD FUNC 1/40 



NCUC(7B-12, 12, 35, 65, x) = 1396.907267 
ALUC(75-12, 12, 30, 35, 65, x) = 69S4.53633S 
1396.907267-5 = W2&WB&H3S 


Example 14: Projected Unit Credit (PUC) 

The PUC cost method adds the use of a salary scale to the TUC method. The PUC assumes that the 
pension benefit earned in any year is based on the final salary or final average salary. Final salary is 
the salary in the year preceding retirement. 


X age 40 earns a pension benefit each year equal to 2% of his final salary. His current salary is 
$45,000 and is expected to increase by 2.5% a year. He was hired at age 30 and expects to retire at 
age 65. What is the NC of his benefit earned in the following year and what is the AL? 


The NCUC() and ALUC() functions along with the final average salary function, FAS(cs, g, cage, 
rage, yrs). ‘cs’ is the current salary, ‘g’ is the expected growth rate of the salary, ‘cage’ is the current 
age, ‘rage’ is the retirement age, and ‘yrs’ is the number of year’s salary to be averaged. Use 1 for 
‘yrs’ if only the final salary is to be used. 



Fi7»nriT fet Y f::t y fht y F? t i _ 
t j— |fl 1 gebr a |Ca 1 c |Qiher |Pr gn I Q |C 1 ean Up 


D 


' f as(45QGQ , 2. 5 , 4G , 65 , 1) 81392.667731 

' 31392. 667731217 -2* 1627.353355 

' ncucC 1 627 . 35 , 1 2 , 43 , 65 , >0 3397 . 74942G 

' ncuct^ ■ fas(45GQQ , 2. 5 , 4G , 65 , 1) , 12 , 40* 

3397.756422 


■ BOB. 2 .5 .40.65 .1) .12 .40.65 .x)l 

ACTUARY RAD AUTD FUNC 4/10 


The final salary is $81 392.67. The pension benefit earned in the following year is 2% of the final 
salary or $1 627.85. The NC is therefore $3 397.76. 
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Everything can be put together at once. 
The AL is $33 977.56. 


m 


|fllgebra|Calc|Oiher|PrgpiIO|Clean 




■ f as(450QG , 2. 5 , 40 , 65 , 1) 31392. 667731 

■ SI 392. 667731217 ■ 2^ 1627.353355 

■ ncuc( 1627. 85 , 12 , 4G , 65 , x) 3397. 74942G 

■ ncuc(2^ ■ fas(45QQQ , 2. 5 , 4Q , 65 , 1) , 12 , 4Q* 

3397.756422 

■ aluc(2^ ■ fas(45GGG , 2. 5 , 4G , 65 , 1) , 12 , 30* 

33977. 56422Q 

... ,2.5, 40, 65, 1>, 12, 30, 40, 65, x>l 

ACTUARY RAD AUTD FUNC E/30 


Example 15: Entry Age Normal (EAN) Cost Method 

The Entry Age Normal normal cost is such that at age e, the entry age into a pension, the present value 
of all future normal costs equals the present value of all future benefits. The entry age normal method 
requires an estimate of a person’s full retirement benefit. 

X age 45 plans to retire at age 65 and receive a pension of $500 a month. He was hired at age 35. 
What are the NC and AL? 


To compute the NC use the NCEANLD(br, mth, eage, rate, person) function where eage is the entry 
age into the pension plan. The LD stands for level dollar; i.e. the normal cost will be the same amount 
every year (unless assumptions change and then a revised normal cost is computed). (If assumptions 
change and normal costs are revised, a supplemental liability may be created which may be funded 
through supplemental costs added to the normal costs.) 


The normal cost is $650.57. 


( FiT Y' fT 

[ToolsJjabl 


t V nr 
lesjps & 


^ Y Y feT 

< Qs|flHHUI TIES|flDU 




ALEANLD(br, mth, eage, cage, rage, person) is used 
to compute the actuarial liability. 


The AL is $9,158.44. 


■ nceanld(5QQ- 12, 12, 35, 65, x) 65G.57 

■ aleanldCSGG- 12, 12 , 35 , 45 , 65 , x) 

9158.44 

aleanld<5QQ*12 ,12 ,35,45,65 ,x)l 

ACTUARY RAD AUTD FUNC E/30 


The Entry Age Normal cost method usually uses a salary increase function when the benefit is based 
on a final salary or final average salary. When salaries are assumed to increase, the normal cost is 
defined as a level percentage (LP) of salary. When this is the case NCEANLP(br, mth, sc, eage, cage, 
rage, person) and ALEANLP(br, mth, sc, eage, cage, rage, person) are used to calculate the NC and 
AL. 


X age 35 plans to retire at age 65 and receive 1% of his final monthly pay rate per year of service. His 
current salary is $2,000 per month and he was hired at age 25. His salary increases 5% each year. 
What is the normal cost and actuarial liability? 

X’s final average salary is $98,787.25, 

His annual retirement benefit is $39,514.90 (1% of final salary times forth years of service). 

And the NC for the year is $1,700.50. 

This represents 7.085 percent of his current salary. Assuming salary projections do not change, the 
NC of future years will be 7.085 percent of the projected salary for that year. 

The AL is $18,008.99. 


You can check all calculations on the next screen shots. Needs some calculation time on the handheld. 
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FAS(2000-12, 5, 35, 65, 1) = 98787.25428 
98787 . 25428 -IS- 40 = 39514.90171 
NCEANLP(39514 ,9 , 12, 5, 25, 35, 65) = 1700.500461 
1700.500461 

= 0.07085418587 

2000-12 

ALEANLP(39514.9, 12, 5, 25, 35, 65) = 18008.98627 


r fit y FT? y fT? Y pT? t fT? 
|Tools|Tables|Fs k Qs|flHHUITIES|flDU RHHh | | 


■ fas(20Q0 ■ 12 , 5 , 35 , 65 , 1) 98787.25 

■ 98787.254289158 ■ l* ■ 4Q 39514. 9Q 

■ ncean 1 p( 395 14.9 , 12 , 5 , 25 , 35 , 65 , x) 

17QQ.5Q 

_ 17QQ. 5QQ4613301 ___ 

2QQQ12 

■ aleanlp(39514. 9 , 12 , 5 , 25 , 35 , 65 , x) 

18QQ8. 9883 


...nlp<39514-9 ,12,5,25,35,65 ,x)l 

liLTUliFiV FiliD HUTU FUhK 5/30 


3. Life Insurance 

Premiums for life insurance may be calculated with 
the functions under F6 and then F3: LIFE INS. 


The basic life insurance functions are Cx(age, person), Mx(age, person) and Rx(age, person). Mx() is 
the sum of Cx()’s from age to the end of the mortality table. And Rx() is the sum of Mx()’s from age 
to the end of the mortality table. 

Cx() is defined as (l x - l x +i)*v x+1 , where l x is the number of persons living at age x. Cx() therefore 
represents the number of persons dying from age x to age x + 1 times an interest rate factor evaluated 
at age x + 1 . It is used to calculate the value of life insurance paid at the end of the year of death. Cx() 
is rarely used in life insurance calculations because most people buy more than just a one-year term 
insurance. Mx() is mostly used. 

Example 16: Whole Life Insurance 

What is the APV or net single premium of a whole life insurance policy of $100,000 payable at the 
end of the year of death for Y age 25? The APV of a whole life policy is given by Mx/Dx. 

The APV of the policy is $4,460.37. 

If Y wanted to pay annual premiums on this policy at the beginning of each year for as long as she 
lives, what would they be? 

The annual premiums would be $264.26. Notice that Mx is divided by Nx where Nx represents a life 
annuity. 

Monthly premiums payable at the beginning of each period would be $22.64. 

100000- M(25, y) 

= 4460.370167 

0(25, y) 

100000- M (25, y) 

= 264.2607926 

N(25, y) 

100000- M (25, y) 

= 22.64336979 

12- NM(25 , 12, y) 


r FIT Y F£t Y F:=t 

|Tools|Tables|Ps & 

QsIfiNNUI T I EsInDU 5 FINnI^TI 

1QQQQQ ■ nx( 25 , y) 
dx(25 , y) 

B 1QQQQQ ■ nx( 25 , y) 
nx(25 ? y) 
1QQQQQ ■ rox(25 , y) 

446Q.37Q2 
264. 2603 
22.6434 

12-nxm(25, 12, y) 

... 0#inx<25, y>/<12#nxm<25,12 

ACTUARY FiliD liUTD 

FUNC 3/30 


ifl'T F£t I 

MB 


'i 

h Ipensions 

Bu 


1 






3= 1 

Mx< 



4s| 

hixnit 



5a| 

Rx< 



6: 1 

Rxrot 


■ ClrHone 



Done 


1 TYPE DR USE + [ENTER]=DK liND [ESC]=CANCEL S 
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Example 17: Term Insurance 


Y want to insure her life for $1,000,000 during her working years from age 25 to age 65, payable at 
the moment of death. But she only wants to pay annual premiums for 20 years. What is the amount of 
the annual premium? For this is example use the Mxm(age, mth, person) function with mth set to oo. 
For life insurance, mth determines when the benefit is paid. A positive number means that it is paid at 
the end of the period of death. A negative number means that it is paid at the beginning of the period 
of death. For example, if mth = 12, the death benefit is paid at the end of the month of death. If mth = 
1, it is paid at the end of the year of death. And if mth = oo, it is assumed to be paid at the moment of 
death. If you noticed that the sign of mth means the reverse for life annuities, you are correct. Again, 
this is a consequence of the mathematical definition of the actuarial functions. 

( fit y FT^ Y fyY y pF y fFP '|f?n 

|Tools|Tables|Ps fc Qs|RNNUITIES|RDU RHHk I I 


1000000- (MM (25, oo, y) - MM (65, oo, y)) 
N(2S, y) - N(45, y) 


1230.800745 


The annual premium is $1,230.80. 


1QQQQQQ (nxK25 ? , y)) 

" nx(25 ? y) - nx(45 ? y) 


1230. S007 


■■■65 J .M J .y>>/<nx<25 J .y>-nx<45 J .y>> 

ACTUARY F;mD hUTD FUNC 1/30 


Example 18: Increasing Insurance 

What is the APV of a life insurance policy that pays $10,000 the first year if death occurs and in- 
creases by $10,000 each year after to X who is age 35? What about Y who is age 35 (DNL#57)? 


10000- R(35) 

= 33520.07425 

0(35) 

10000- R (35, y) 

= 28970. 1B447 

0(35 , y) 

The APVs are $33 520.07 and $28970.10. 

5. Non-integral Ages and the Uniform Distribution of Deaths 

All the commutation functions accept non-integral ages, e.g. 50.5 or 32.3. Their values are calculated 
on the assumption that deaths during the year follow a uniform distribution (the variable Extrapol() is 
set to udd). However, if you want a straight liner interpolation between non-integral ages, set Ex- 
trapol() to lin. If non-integral ages are used with function such as Nxm() or Mxm() where 
abs(mth) > 1 and assuming udd, there are small relative errors in the calculations. Towards the end of 
the mortality table the errors can become large. If non-integral ages are used only when abs(mth) = 1 
or if only integral ages are used when abs(mth) > 1, there is no error given the uniform distribution of 
deaths assumption. 


I^TpENS FoNsTl I FE ’i NsT 

1 

b 10000-rx(35,x) 
dx£35 , x) 

33520.0743 

1000G ■ r x( 35 ,y) 
dxC35 ? y) 

2S970. 1045 

10000#Rx<35 _y>/dx<35 
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FUNC E/30 
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7. Q’s and P’s 

The Q’s and P’s associated with actuarial math are the 
probabilities of dying and living. (In actuarial textbooks 
they are always lower case. I use upper case because the 
TI-89 does not do subscripts.) The Q’s and P’s menu is 
under F3. 


lx(age, person) returns the number of persons living at the beginning of age x for person X, Y or J. 
The ages extend from age a to age cd, the beginning and ending ages of the current mortality table. 
The number of persons at age a is the number chosen for the radix of the mortality table, in this case 
1,000,000,000. For instance, the number of persons X living at the beginning of age 65 is 
866,252,268.952. 

Px(age, person) gives the probability of a person of age x living one year. For Y age 25 the probabil- 
ity of living one year is 0.999687. Px() = l x +i/l x 

Qx(age, person) gives the probability of a person of age x dying in the next year. For Y age 65 the 
probability of not being alive at age 66 is 0.009286. Qx() = (l x - l x +i)/l x 

crm(age, person) returns the central rate of mortality of a person of age x. For Y age 65 the central 
rate of mortality is 0.009329316014. crm() = (l x - l x +i)/LLx(age, person) where LLx() is defined be- 
low. 

crm for DERIVE is new! 

~lx(6B) = B662B2268 .9 
Px(2B , y) = 0.9996B69999 
Qx (6B , y) = 0. 0092BB999913 
crm(6B , y) = 0. 00932931B926 


r FIT Y Y F5t Y FHt Y F£t ^fiY 'i 

|Tools|Tables|Ps & Qs|flNNUITIES|flDU RNNb | | 

■ lx(65 7 x) 

866252268. 952 

■ px(25 7 y) 

.999637 

j* 

in 

■jj 

X 

V 

■ 

.009286 

■ crn(65 ? y} 

.009329 

crnK65, v> 

ACTUARY RAD liUTD 

FUNC H/30 
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■ ClrHone 

2: Fx(^^ 
3: Qx( 

4: crnt 
5: nPxt 
6: nQxt 
7:tnQxC 
8:LLxC 
9: Tx( 
R:elxC 

[jETOE 

ClrHome 

TYPE DR USE + [ENTER]=DK AND CESC]=CANCEL 


nPx(age, n, person) gives the probability of a person of age x living n years. For X age 34 the prob- 
ability of reaching age 58 is 0.945032212698. 

nPx() = l x+n /l x 

nQx(age, n, person) of course is the probability of a person of age x dying before reaching age x+n. 
For X age 34 the probability of dying before reaching age 58 is 0.054967787302. 

nQx() = (l x - l x+n )/l x 

tnQx(age, t, n, person) returns the probability of a person of age x living t years and then dying in the 
next n years. For Y age 34 the probability of living 10 years and they dying in the next 5 years is 
0.005590571447. 

tnQx(age, t, n person) = nPx(age, t, person) * nQx(age+t, n, person) 

LLx(age, person) gives the number of life-years lived by the l x persons who attain age x over the year 
from age x to x+1. On the assumption that deaths are uniformly distributed, LLx() = (l x + l x +i)/2. For 
example, the number of life-years lived over the year by the X persons who attained age 45 is 
970,159,731.358 years. 
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Tx(age, person) returns the total future lifetime of the l x persons who attain age x. The total future 
lifetime of the X persons who attained age 45 is 3,367,321,573.8 years. 

Tx() = l x /2 + l x+1 + . . . + 1 (0 

elx(age, person) gives the complete expectation of life of a person age x. For Y age 65 the expected 
life is 20.69 years. 

elx() = Tx()/lx() 


nPx(34, 58 - 34) = 0.9450322126 
nQx(34 , 58 - 34) = 0.05496778731 
tnQx(34, 10, 5, y) = 0.005590571399 
nx(45) = 970159731,3 
tx(45) = 33673281571.2 
six (65, y) = 20.69167200 


fFiTurfr FF? Y f:=t t fht y ft Y Fe^ Y 'i 
| t fl 1 gebra |c a 1 c lother |Pr gm 1 0 |c 1 ean Up | | 
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■ crn(65 ? y} 

.009329 

■ npx(34, 53 - 34 , x) 

.945032 

■ nqx(34 ? 53 - 34 ? x} 

. 054968 

■ t-nqx<;34 7 10 ? 5 ? y} 

.005591 

■ 1 lx(45 ? x} 

970159731.35S 

■ t-x(45 ? x) 

33673281573.8 

■ elx(65 ? y} 

20.691672 

elx<65.w>l 

liLTUhRY F;hD hUTD 
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7. Conclusion 

I hope you’ve enjoyed this brief excursion into actuarial mathematics. This is only a taste of what the 
field has to offer. If you have any questions, comments, suggestions please email me at 
don.phillips@gmail.com . 



I have the pleasure to present the program committess for the two strands of this conference. 
Co-Chairs are Michel Beaudin & Josef Bohm / Berhard Kutzler & Vlasta Kokol-Voljc. 
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TI-Nspire CAS = Successor of Derive!? 

by B Kutzler 


TI-Nspire CAS is a powerful new mathematics tool suited for teaching and learning. It combines several appli- 
cations under one roof: A computer algebra system (CAS), a graphing software, a dynamic geometry system, a 
spreadsheet, a text editor, and a data & statistics package. 



The Calculator application is a computer algebra system (CAS) based on the code of the handheld Voyage200. 

The Graphs & Geometry application is a double application. It is a graphing application as known from graphing 
calculators such as the TI-83 or TI-84. It is also a dynamic geometry application as known from dynamic geo- 
metry products such as Cabri Geometry. Having both applications integrated with each other creates a powerful 
“geo-graphical” environment which goes far beyond what the individual applications usually offer. 

The Lists & Spreadsheet application is a spreadsheet as known from programs such as Microsoft Excel com- 
bined with a list editor as known from the Voyage200. 

The Notes application is an easy-to-use word processor which allows modes and styles particularly useful for 
annotating mathematical and educational work. 

The Data & Statistics application is an easy-to-use package for statistical analysis with features as known from 
programs such as Fathom. 

The seamless integration of these key technologies is one of the major strengths of TI-Nspire CAS. The indivi- 
dual technologies have been on the market for quite some years and have been tested in the classroom. Litera- 
ture is available on that and there is strong evidence for each of these technologies that their proper integration 
into teaching and learning is beneficial for students. TI-Nspire CAS integrates all of these technologies into a 
powerful mathematics teaching and learning environment. This integration of well established technologies into 
one tool is a quantum leap and takes us into a new dimension of supporting students in understanding and learn- 
ing mathematics. 

TI-Nspire CAS is the proclaimed successor of Derive. Some Derive users deny this with the argument, that TI- 
Nspire CAS still lacks several of the Derive features. Is TI-Nspire CAS the successor of Derive - or is it not - or 
can it become its successor? I consider TI-Nspire CAS the successor of Derive and will provide my arguments 
in this article. 

Some people agonize about the fact that Derive is not further developed. What does the history of Derive have 
to say about this? 
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Soft Warehouse was founded 1979 by Albert D Rich and David R Stoutemyer as a partnership, later it got 
incorporated as Soft Warehouse, Inc. with offices in Honolulu, Hawaii. Soft Warehouse released their first CAS 
product in 1979. It was called muMATH-79 and ran on 8080 and Z80 computers with as little as 48K bytes of 
memory running the CP/M operating system. Later versions were muMATH-80 and muMATH-83. 

In October 1988 Soft Warehouse, Inc. released Derive as the successor of muMATH. Derive was a DOS 
program running on PC compatible computers with a minimum of 512K bytes of memory. In 1996 Derive for 
Windows was released. It had a GUI Windows interface, a 32-bit math engine kernel, and ran on PC compatible 
computers running MS Windows and NT. In 1999 David and Albert sold their company, and with it Derive, to 
Texas Instruments Incorporated. 

The authors of Derive wrote in their “ Brief History of the muMATH/Derive CASs” the following sentence: 
“Rather than just refining and improving muMATH, we decided that an entire re-write was needed. Derive is 
the result. Instead of being written in muSIMP, Derive is written directly in LISP, specifically muLISP. ” 

After 9 years of muMATH the authors decided that it is time to do a quantum leap into a new technology by 
rewriting their CAS know how. 

After 19 years of Derive it seems reasonable to do the same. Derive has been a key to changing the world of 
teaching and learning mathematics. Derive has changed the life of many teachers and students and has triggered 
many new developments, including the development of the legendary TI-92 handheld. But, the design of Derive 
is from the 80s, which makes it more and more difficult to keep up with the more and more sophisticated desires 
of users in the new millenium. It is time for a quantum leap into a new technology by doing another rewrite. And 
this has happened with the development and release of TI-Nspire CAS. 

Some people agonize about the fact that some features of Derive are not (yet) available in TI-Nspire CAS. 

muMATH had some nice features which were not available in Derive. For example, muMATH was an open 
source code product written in muSIMP. The user could look at the programs and functions and could write their 
own programs and functions. Derive was written in muLISP . The program code was not available any more to 
users, nor could they do any programming in the first version(s) of Derive. 

It is true, that several Derive features are not (yet) available in TI-Nspire CAS. But it is also true that there 
are very many powerful features in TI-Nspire CAS which were not available in Derive, although quite a few of 
them were on the “Derive Wish List” for many years - and probably some of these features we never would 
have seen in future versions of Derive, simply because the design of Derive did not allow for their implementa- 
tion. Later in this article we will look at the “Derive Wish List” and look at features on this list which are al- 
ready available in TI-Nspire CAS. 

In April 2007 a “Derive-TI-Nspire Transition Conference” was held in Austria with about 20 Derive experts 
from all around the world working with about 10 people from the TI-Nspire CAS development team, including 
David Stoutemyer. The goal of this meeting was to pin down those Derive features which should go into future 
versions of TI-Nspire CAS. 

In 1988 

- Derive was the proclaimed successor of muMATH, 

- Derive was an entire re-write in another language, 

- Some features of muMATH were not available in Derive Version 1. 

In 2007 

- TI-Nspire CAS is the proclaimed successor of Derive, 

- TI-Nspire CAS is an entire re-write in another language, 

- Some features of Derive are not available in TI-Nspire CAS Version 1. 

History repeats! 

A famous quote says: “The revolution eats its children ”. This has proven true throughout history - and mathe- 
matics is not any different from other parts of life (even if non-mathematicians often doubt that; but, yes, mathe- 
maticians are humans, too!). While we see changes of products, there is also continuity. The continuity facor 
here is David Stoutemyer. David is the author of muMATH, David is the original author of Derive (until Albert 
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Rich took over, who is the original author of muLISP, the language in which Derive was written), and David is 
the original main author of the CAS code of TI-Nspire CAS. So “David’s CAS Code” is the “immortal soul” 
which already sees its third reincarnation. Isn’t this a proof that reincarnation exists? Why grief about death 
when evidently there is resurrection? The new “body” looks different, the new “body” speaks a different lan- 
guage, the new “body” sometimes acts different, but the new “body” has the same soul! The new body is very 
young and very strong, lives in a promising and supportive environment, and definitely has got everything it 
needs so that it can learn a lot and grow stronger and stronger in the years to come. 

In the previous 19 years the activities around Derive created the “Derive Spirit”. This spirit lives through people, 
hence it is still alive, it will stay alive, and it will continue to improve the way we teach and learn mathematics. 
The embodiment of the “Derive Spirit”, the “Derive community” is united by the vision to use CAS for making 
mathematics more accessible to students. This vision is a pedagogical version of the original vision of the Derive 
authors “ to make computer algebra widely available to the masses on small computers 

One of the Derive community’s main propositions was that the tool is only a means to reach the goal. We 
always advocated the use of Derive when its use supports the reaching of a (given) teaching goal - one should 
not use Derive (or any other tool) for its own sake. Let’s stick to our vision - rather than to a product! 

A glass half filled with beer can be considered half empty or half full, dependent on whether your attention is 
on the full half or on the empty half. It’s the same here. The change from Derive to TI-Nspire CAS can be seen 
as having lost Derive or as having won TI-Nspire CAS - with all its current capabilities and with its potential for 
the future. We go out to teachers and preach innovation. We ask them to allow for changes in their way of teach- 
ing. Now we are challenged to allow for a change of the product, which, in fact, is a change to a new format of 
“David’s CAS Code”. 

The Derive community has a great history of influencing (could we say “d(e)riving”?) the further develop- 
ment of Derive. Having had to speak to a small team of only two programmers made communication easy and 
decision chains short. Sometimes we saw a bug fix or a new feature over night - while other features never 
made it, because we were not “convincing enough” . . . Small teams have their advantages - and large teams 
have their advantages too. Clearly, in a large team decisions take longer, but with a large team there is more 
“power” behind development, so that more wishes can be fulfilled. And there is also more “power” behind mar- 
keting the product, including hot line support, teacher training, etc. TI-Nspire CAS definitely can become much 
bigger than Derive ever could have become. And this is good for our vision . . . 

During my years of being responsible for Derive, I compiled a list of features which Derive users from around 
the world had asked for. Most items on this list are for the user interface and the interface to other technologies. 
Relatively few items are for changes to the CAS engine. This reflects the fact that many users thought that the 
CAS engine was pretty satisfactory for their purpose, they just wished for more comfort and more flexibility in 
connecting to other tools. In fact, quite a few items from this Derive Wish List would have been very difficult, if 
not impossible to be implemented in Derive. Some items from this list are already available in TI-Nspire CAS, 
others could well be seen in future versions, because the new design of TI-Nspire CAS allows for these features. 
Following is a selection from the Derive Wish List and my comments regarding their availability in TI-Nspire 
CAS. 

“ Manipulate the graphics window with the mouse (change scales , move origin , . . .) 

The TI-Nspire Graphing application offers a lot of options for changing the window settings. 

- You can grab the plane and shift it, thus panning it within the window: 
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-You can click on an axis’ maximum or minimum value, and then change the value via the keyboard: 


X 

X 

X 



12 


-You can use the Window Settings command from the Window menu, and then enter the values you want. 




Window Settings 



Zoom - Box 


P 

Zoom - In 


P. 

Zoom - Out 



Zoom - Standard 



Zoom - Quadrant 1 


l-ft 

Zoom - User 



Zoom - Trig 



Zoom - Data 



Zoom - Fit 


-You can use one of the many zoom commands from the Window menu. 


-You can grab any of the tick marks of an axis, and then change the scale factor by dragging it: 
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-You can click on the axis, and then right click for the context sensitive menu. The Attributes command 


opens a graphical menu in which you can change the appearance of the axis (display it with or without ar- 
rows). This menu also offers previews to the most important zoom commands from the Window menu. 


^ Recent ► 
Attributes 
Hide Axes 
Show Grid 
Zoom ► 


2 

' + 

% 



(3/3) No arrows 
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“ Data exchange with Cabri products ”: 

Cabri Geometry is the technological basis of the Graphs&Geometry application, so basically now you have 
both technologies (and more) within one product. 

“ Simple drawing options with ruler , compass , and protractor”: 

This is exactly what the Graphs&Geometry application offers through the Cabri technology. 

“ More links between the algebra and the graphing window (such as transferring graphical data)”: 
In an TI-Nspire Graphing or Geometry window you can store a value such as a coordinate or a measurement 
(of a length, an area, a slope, an angle, or an integral) in a variable. 
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Store Var 

Unlink 
LinkTo: ► 


(1 | := - 9.18681 

, 4 . 48352 ) 

• 




This variable and, hence, its value, can be used in any other TI-Nspire application. 

Alternatively, you can link a value in a Graphs&Geometry text box to any variable which was defined in any 
other TI-Nspire application. 


Store Var 

Unlink 


LinkTo: 




W 



This variable/value can be transferred to geometrical objects as a distance. Hence, distances in a Graphing or 
Geometry window can be controlled from any other application. 


“Update the algebra expression when moving the graph”: 

In TI-Nspire you can grab a curve and change it. The function’s equation will change with it. 
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This feature is available for the following classes of functions: 

- linear function y-b 

- linear function y = a • x + b 

- quadratic function y - a • (x - b ) 2 + c 

- exponential function y = e (a ' x+b) + c 

- exponential function y = b • e ax + c 

- exponential function y = d- e a ' x+b + c 

- logarithmic function y = a- In (c • x + b) + d 

- sinusoidal function y = a- sin(c • x + b) + d 

- cosinusoidal function y = a- cos(c • x + b) + d 


“Animation of graphics”: 

The Graphs&Geometry application offers a tool for animation. This tool can be applied to “independent” 
points on an object such as a segment, a ray, a line, a circle, or a curve. By the nature of the dynamic geome- 
try application, all objects dependent on this point will be animated together with the point. 





A related tool is Locus. It draws a family of objects generated by the movement of an object it depends on. 
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“Dynamic interactive graphics (change a point and the regression line changes )”: 

In the Data&Statistics applications you can grab points of a data plot or a scatter plot, move them, and then 
see the changes of dependent objects such as regression lines or regression curves. Also the numeric repre- 
sentations, for example in the spreadsheet, change accordingly. 



“Import Excel data”: 

In the Lists&Spreadsheet application you can paste data which previously was copied into the clipboard from 
an Excel document. 
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“Dynamic algebra documents (also called 6 history refresh when you change an algebraic expression, then 
all objects depending on it shall be updated) 

The TI-Nspire Calculator does not offer such a feature, but one can build a problem solving template in the 
Lists&Spreadsheet application (in which most Calculator/CAS features are available). Therefore by the na- 
ture of a(n algebraic) spreadsheet, this gives a dynamic algebra document. The following two screen images 
show a template for computing the coordinates of the point of intersection of two lines given by their equa- 
tions. 


1 first line 

2* x +3*y=5 


\2 second line 

8*x-13*y=21 


3 first variable 

X 


4 second variable 

y 


5 

x=64/25 and y=... 


6 x-coord intPnt 

64/25 


7 y-coord intPnt 

-1/25 



1 first line 

2*a*x+b*y=? 


2 second line 

b*x-3*a*y = 2 


3 first variable 

x 


4 second variable 

y 


5 

x=(21 *a+2*b)/(6.„ 


6 x-coord intPnt 

(21 *a+2*b)/(6*a... 


|? y-coord intPnt 

-(4*a-7*b)/(6*a... 



The current version of TI-Nspire CAS still has limitations which prohibit an extensive exploitation of this 
very powerful capability. For example, cells cannot contain lists or matrices. But the mechanism is there, so 
all it needs is removing some of these limitations. 


Other wishes which are already fulfilled in the current version of TI-Nspire CAS include: 
- “Compute all solutions of trigonometric equations ” 
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- “Offer a simple function for solving ordinary differential equations 
de So 1 ve (y =>+ x, x,y) y=c1-e x -x~l 


It always takes courage to do quantum leaps such as from muMATH to Derive or from Derive to TI-Nspire. I 
guess there is nobody in the “audience” of this article who would oppose the statement that the quantum leap 
from muMATH to Derive was a great step forward and definitely worth having been done. This quantum leap 
also resulted in a significant increase of users. There are many more Derive users than there were muMATH 
users. With the technical and commercial power of Texas Instruments behind TI-Nspire CAS it is very likely 
that soon there will be many more TI-Nspire CAS users than there are Derive users. And isn’t this what we 
wanted to achieve with all our activities of promoting the use of CAS for teaching and learning mathematics? 

In a follow-up article I will demonstrate how easy it is in TI-Nspire CAS to combine its applications to create a 
“holistic” experience of solving an optimization problem geometrically, graphically, and algebraically. 


As Bernhard mentions in his last paragraph he sent a second paper demonstrating the 
Nspire. At this occasion I'd like to invite all of you who have been working and experimenting 
with the Nspire (either handheld or PC-version) to share your findings with us. I also repeat 
asking for TI-92+, TI-89 and V 200 papers. Just recently checking Nils Hah nf eld's CME (see 
page 3 and reproducing Josef Lechner's “Falling Bodies” I was again fascinated by the 
power of the handheld. Come on and share your treasuries (and your problems, of course). 

I send my “hand held regards”, 

Josef 


I met David Sjostrand in March when we discussed which features of DERIVE we 
would like to find in future versions of TI-Nspire. One of many points was, that we 
miss ITERATE and ITERATES as a powerful tool - even in times of program- 
ming, because of its didactical value. David presented a fine example and I 
asked him leaving it for the DNL. He agreed without hesitating and he sent an- 
other paper dealing also with triangles: 


Hi Josef, 

here is my paper again. I would be interested to know if Corollary 1 is already known. I didn't know 
this until I was working with this paper. Perhaps you or some of the readers of the DNL know. 

Best regards, 

David 
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Incenter, Excenters, Orthocenter, Centroid and 
common Intersection Points in a Triangle 

David Sjostrand, Onsala, Sweden 

Let ABC be a triangle. In this paper we often identify a point P with the vector OP . 

Definition 1: 

1 ine(A, B , x , y) := (A -B)-(x-A)-(A - B ) ■ (y - A ) = 0 

2 2 111 2 

D is a point lying between the points A and B on the line passing A and B dividing the segment AB in 
two parts in the ratio alb , where a,b> 0 , counted from B. 

7 / __ x / __ . \ aA + bB aA + bB 

Then b(B -D) = a(D- A) D = . It is obvious that the point is on the line 

passing A and B for any real numbers a and b such that a + b + 0 

Definition 2: 

A and B are two points in the plane and a and b are real numbers a and b such that a + b + 0 , then the 
aA + bB 

point D = divides the line segment AB in the ratio alb counted from B. 

a + b 

In Defintion 2 D is an arbitrary point on the line passing A and B. 

We will use Derive to prove 


Theorem 1 

The lines 

( bB + cC \ ( cC + aA \ ( aA + bB \ 

Line A , ,x,y , line B , ,x,y and line C, ,x,y are concurrent and 

V b+c ) V c+a ) I a+b ) 


their common point of intersection is 


aA + bB + cC 
a + b + c 


Proof: 

We make the following definitions in Derive 

[A == [al, a2], B == [bl, b2], C == [cl, c2], a == a, b == b, c == c] 

If we simplify 


T f b-B + a-A A f a-A+c-C 

SOLUTIONS line C, , x, y , line B, , x, y , [x , y] 

b + a J L a+c JJ 


SOLUTIONS line A 


^ r a- A + b-B \] 

x, y , line C, , x, y , [x, y] 

J L a+b JJ 


we receive in both cases 


a-al + b-bl + ocl a-a2 + b-b2 + c-c2 


a + □ + c 


a + d + c 


which is equal to 


aA + bB + cC 
a + b + c 


Q.E.D. 
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There is a converse of Theorem 1 : 


Theorem 2 

A x is on the line passing B and C, B x is on the line 
passing C and A and C x is on the line passing A 
and B. 

If the three lines line ( A, A x ), line (B,B X ^ and 
line(C, C x ) are concurrent, there are three real 

numbers a, b and c, such that 

, bB + cC _ cC + aA aA + bB b 

A x = , B x = and C, = 

b+c c+a a+b 



Proof: 


Figure 1 


Since A x is on the line passing B and C there exist real numbers b and c such that A x = 


bB + cC 


Since B x is on the line passing C and A there exist real numbers c x and a x such that 


— (c x C + a x A) 

„ c,C + a,A c, cC + aA c 

B x = — — = — = where a = — a x . 

c ' +a . -h+o.) c+a c i 

c v 7 

C 1 

( aA + bB\ 

The lines lineiA, A ), line(B,B x ) and line C, are concurrent according to Theorem 1. 

I a + b 


Therefore C = 


aA + bB 


Q.E.D. 


Corollary 1 

AC BA CB 

The lines line (A, A 1 ) , line^B, B x ) and line(C , C, ) are concurrent if and only if = 1 

Cj B A x C B x A 

Proof: 

It follows from Definition 2 and Theorem 1 and Theorem 2 that lines line(A,A x ^, line(B,B x ^ and 

line(C , C, ) are concurrent if and only if there are real numbers a , b and c such that ^ + divides 

BC in the ratio b!c counted from C. 


cC + aA 

divides CA in the ratio c! a counted from A. 


c + a 
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aA + bB 


divides CA in the ratio alb counted from B. 


n . AC, BA, CB, b c a 

Since L 5 L = = 1 Corollary 1 is proved. 

CB AC B,A cab 


Corollary 2 

If a = \B - Cl , b = \C - A\ and c = \A — B\, the incenter of the triangle ABC is 


Q.E.D. 


aA + bB + cC 
a + b + c 


Proof: 

bB + cC 


- divides BC in the ratio b/c counted from C. 


cC + aA 


divides CA in the ratio c/a counted from A. 


aA + bB 


divides AB in the ratio alb counted from B. 


( bB + cC \ 

If a = |Z?-C|, b = |C — ., 4 | and c = \A -B\ then the lines line A,— ,x,y , 

^ cC “l - ^ f ayi “I - bB 

line B, , x, y and line C, , x, y are the angle bisectors of the triangle ABC. 

y c + a J y a + b J 

The incenter of a triangle is the common intersection point of the angle bisectors of the triangle. 

Now the Corollary 1 follows from Theorem 1 . 

Q.E.D. 


Lemma 

Let ABC be a triangle and 
a = |i?-C|, & = |C — ^4| and c = \A — B\ 

If AD is the bisector of the exterior angle 
„ —bB + cC 


BAF, then D = ■ 


-b + c 



Proof: 

BE is parallel to CF. Therefore BE = BA and the triangles BED and DCA are similar which gives 

CD b J 

= — and consequently: 

BD c 


'.(D-C) = b(D-B)<*D = 


-bB + cC 
-b + c 


Q.E.D. 
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Corollary 3 

The excenters of the triangle ABC are 

-aA -bB + cC aA-bB-cC 

, and 

-a-b+c a-b-c 

-aA + bB-cC 

-a + b-c 



Proof: x 

f -fog + C c \ 

The excenter, F, is the intersection point of the bisector of ZDBA , line A, , x, y and the 

V ~b + c J 

^ £ 7/4 + cC ^ 

bisector of Z.BAF , line B , ,x,y . 

{ -a + c ) 


7 . ( A -bB + cC ) 

According to Theorem 1 the intersection point of line A , , x, y and 

V -b + c J 

f -aA + cC V -aA-bB + cC 

line B , ,x 9 y is . 

I -a + c ) -a-b + c 


We have to take care of the case that b = c. In this case the bisector of the exterior angles of A is paral- 
lel to BC. The equation of this bisector is 


(C - B).0-A)-(C - B )'(y-A) = 0 

2 2 111 2 

You receive the excenter as the intersection point between this bisector and 

f bB + aA \ . 

line C , , x, y using Derive : 

V b + a ) 


SOLVE line C 


a-A + b-B 


y , (C - B - A) - (C - B My - A ) = 0 , [x , y] 
J 2 2 1 11 2 


i-al + b-(bl - cl) a-a2 + b-(b2 - c2) 

a y = 

a a 


a- a\ + b(b\-c\^ a- a2 + b(b2-c2) a-al + b-bl-c-cl a-a2 + b-b2-c-c2 
a 9 a L a+b-c ’ a+b-c 

-aA -bB + cC 
-a-b+c 

By symmetry we conclude that the formulas for the excenters are as we have stated in Corollary 2. 

Q.E.D. 
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Using ITERATES for finding the midpoint of a Triangle 

Plot a sequence of triangles, Tn, where the vertices of Tn are the midpoints 
of the sides of Tn-1. TO = P is a given triangle. 

" -1 3 

2 3 

#1: P := 

-2 -2 
_ -1 3 . 

If we make the below defintion of f we have that Tn = f(Tn-l). 

" P + P P+P P+P P + P ' 

12 23 3112 

#2: f (p) := , , , 

L 2 2 2 2 J 

If we plot the below expression we receive the triangles TO, Tl, . . . T20. 

#3: ITERATES(f (x) , x, P, 20) 

The below expression gives the median from A to BC in the triangle ABC. 






D 


-N-L#67 


David Sjostrand: Points of Interest in a Triangle 


p45 


#4: 


MEDIAN (A, B, C) := y = 


f B + C T 


K 2 J 

2 



f BtCl 


k 2 J 

1 



■ Cx - A ) 
1 


+ A 

2 


The below expression gives the medians of TO. 

medians := |"MEDIAN(P , P , P ), MEDIAN (P , P , P ), MEDIAN (P , P , P )1 
#5: L 123 231 3 1 2 J 

The medians of TO are obviously the medians of Tn for all n. 

Let Sn be point in the interior or on the boundary of Tn. 

It is obvius that the sequence Sn is convergent. Let lim Sn = M. It is obvious 
that M is the unique point belonging to all Tn. Since all the medians intersect 
Tn for all n they must all pass M. 

(See the plot of #5 on the next page.) 


Thus we have proved that the medians of a triangle intersect in one point, H. 



Let Pn,l, Pn,2 and Pn,3 be the vertices of Tn. We have that 

Pn , 1+Pn , 2+Pn , 3 = (Pn-1 , 1+Pn-l , 2)/2 + (Pn-1 , 2+Pn-l , 3)/2 + (Pn-1 , 3+Pn-l , l)/2 = 

= Pn-1 ,1+Pn-l, 2+Pn-l, 3 = PO , 1+PO , 2+PO , B . 


On the other hand lim Pn,k = M. 

Therefore 3M = 1 im(Pn , 1+Pn , 2+Pn , 3) = 1 im(P0 , 1+PO , 2+PO , 3) = PO , 1+PO , 2+PO , 3 . 
We have thus proved the wellknown formula M=(A+B+C)/3. 

Let M be the intersection point of the medians. 

We can see that lim Tn = M when n tends to infinity. 

P + P + P 
12 3 

# 6 : 

3 



We can verify that M is given by #6 by verifying a vertice of T20 with #6 


#7: (ITERATE (f(x) ■ x, P, 20)) 

1 


349523 139S103 " 

1048576 ' 1048576 . 


P + P + P 
12 3 

#8: (ITERATE (f(x) , x, P, 20)) 

1 


7 


5 


3 


3145728 


3145728 
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Titbits 34 

or Polynomial Arithmetic for the Advanced 

(c) Johann Wiesenbauer, Vienna University of Technology 

On quite a few occasions I have been dealing with polynomial arithmetic in the past (for 
example, cf. my Titbits 13, 17, 25, 26). As this is a very important and fundamental 
topic, indeed, hopefully you won't object to my resuming it once more. To be more pre- 
cise, I would like to put together a number of basic polynomial routines that could be 
used in a Derive-implementation of the celebrated Schoof-algorithm for point counting 
on elliptic curves of the form y~2 = x~3 +ax +b mod p for some prime p, although this 
algorithm (and some refinements like the SEA-algorithm) are covered in depth only in a 
future issue of this series. Hence, to some extent this is also a continuation of my 
Titbits 30, where I already listed a number of important tools regarding computations 
on those elliptic curves. 

Let's start with the 4 fundamental routines polyrem(u,v,p,x), polyquot(u,v,p,x), 
polygcd(u,v,p,x) and polypower(u,n,v,p,x). In all cases u and v are two polynomials in 
Zp[x], i.e. the polynomial ring over the residue class ring Zp for the prime p. While the 
first two are supposed to compute the remainder and quotient of the division u:v, the 
output of the third and fourth should be the gcd(u,v) and the polynomial power u~k mod 
v, both again in the ring Zp[x], As for the last routine we assume w.I.o.g. that v is 
normed, i.e. has got 1 as the leading coefficient. 

Using basically the routines from the Titbits 25 quoted above, a first implementation 
could look like this: 

poly reni(u , v, p, x) := 

Pnog 

#1: u := FACTOR (REMAINDER (u, v, x), Trivial, x) 

P0LY_M0D(NUMERAT0R (u) ■ I NVERSE_MOD( DENOMINATOR (u) , p) , p) 

po1yquot(u , v, p, x) := 

Pnog 

#2: u := FACTOR (QUOTIENTS, v, x), Trivial, x) 

PQLY_MOD(NUMERATOR (u) ■ I NVERSE_MOD( DENOMINATOR (u) , p) , p) 

po1ygcd(u , v, p, x, r_) := 

Pnog 

u := P0LY_M0D(u , p) 
v := POLY_M0D(v , p) 

Loop 

#3: If v = 0 

RETURN u 

r_ := polyrem(u, v, p, x) 

WRITE (P0LY_DEGREE(r_, x)) 
u := v 
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po1ypower(u, k, v, p, x, u_ := 1) := 

Loop 

If k = 0 
RETURN u_ 

#4: If ODD? (k) 

u_ POLY_MOD(REMAINDER(u- u_, v, x) , p) 
u := POLYJflODCREHAINDERCu-u, v, x) , p) 
k := FL00R(k, 2) 

In fact, they do work, as the examples below show. 

#5: w := P0LY_M0D(3-RAND0M_P0LY(x, 3, 7), 7) 

3 2 

#6: w := 3-x + 6-x + 6-x + B 

#7: u := P0LYJ0D(4-w-RAND0M_P0LY(x l 7, 7), 7) 

10 9 S 6 B 43 2 

#8: u := B-x + x + 3-x + 2-x + 2-x + 4-x + x + 3-x + 6-x + 6 

#9: v := P0LYJ0D(2-w-RAND0M_P0LY(x l 7, 7), 7) 

10 S 7 B 4 3 2 

#10: v := 6-x + 6-x + 3-x + 6-x + 4-x + 3-x + 2-x + 6-x 

3 2 

#U_: polygcd(u, v, 7, x) = 6-x + B-x + B-x + 3 

#12: polygcd(po~lyquot(u , w, 7, x) , polyquot(v, w, 7, x) , 7, x) = 2 

#13: v := POLY_MQD(INVERSE_M0D(6 , 7)-v, 7) 

10 S 7 B 4 3 2 

#14: v:=x + x + 4-x + x + 3-x + 4-x + B-x + x 

98 76 54 3 2 

#15 : po1ypower(u, 10, v, 7, x) = 2-x + 2-x + 5-x + 4-x + 5-x + x + 2-x + 6-x + 6-x + 1 

10 98 76 54 3 2 

#16: P0LY_MOD(REMAINDER(u , v, x) , 7) = 2-x + 2-x + 5-x + 4-x + 5-x +x + 2-x + 6-x + 6-x + 1 

There are a number of objections though. A minor one is that the outcome of 
polygcd(u,v,p,x) is usually not normed, i.e. it hasn't got 1 as the leading coefficient, as 
this is a common convention for gcd's of polynomials. There is an easy remedy for this 
problem, namely by "norming" the divisor polynomial before each division in the euclid- 
ean algorithm, as this was already carried out for v above before the computation of 
polypower(u,10,v,7,x) 

po1ygcd(u, v, p, x, d_, r_) := 

Pnog 

u := P0LY_M0D(u , p) 

If u = 0 

RETURN P0LY_M0D(v, p) 

Loop 

v := P0LY_M0D(v, p) 

If v = 0 

#17 : RETURN u 

d_ := POLY_DEGREE(v , x) 

WRITE (dj 

r_ := LIM(v/x A d_, x, 

v := P0LY_M0D(INVERSE_M0D(r_, p)-v, p) 
r_ := REMAINDER (u, v, x) 
u := v 
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3 2 

#1£ : polygcd(u, v, 7, x) = x + 2-x + 2-x + 4 

A far more serious problem with respect to our goal is the poor performance of the 
routines above, if u and v are polynomials of a "large" degree, say of degree about 1000 
or even more. (For example, the corresponding computations fur polygcd(u,v,p,x) would 
take some quarters of an hour, rather than some seconds, as the interested reader can 
easily check himself by following the pattern above!) 

In the following I'm trying to fix this problem by setting up my own polynomial arithme- 
tic based on the the vector representations of the involved polynomials. In order to 
ease the change between the polynomial representation and the vector representation 
of a polynomial I will first provide the routine polytovec(u,p,x) as well as its inverse 
vectopoly(u,x). 

po1ytovec(u i p, x) := 

Pnog 

u := POLY_M0D(u, p) 

If u = 0 
RETURN [] 

#19 : u := MOD(SUBST(TERMS(u + p- (x A POLY_DEGREE(u) - l)/(x - 1)), x, 1), p) 

Loop 

If FIRST(u) > 0 
RETURN u 
u := REST(u) 

j- 

#20: vectopoly (u , x) := u-VECT0R(x , j_, DIM(u) - 1, 0, -1) 


Furthermore, we need some basic operations for polynomials represented as vectors, 
such as the product vecmult(u,v,p) of two polynomials u and v and the remainder 
vecrem(u,v,p) of the division u:v (as always in Zp[x] !) 


vecmult(u, v, p, s_) := 

Pnog 

If DIM(u) < DIM(v) 

[s_ := u, u := v, v := s_] 

5 := 0- u 

v := REVERSE (y) 

Loop 

If v = [] exit 

#21: s_ := AD]0IN(0, s_ + FIRST(v)-u) 

u := APPENDCu, [0]) 
v := REST(v) 

:= HODC's , p) 

Loop 

:= RESTCO 
If FIRSTCO t 0 
RETURN s_ 


vecrefn(u, v, p, k_ := 0) : = 

Pnog 

If DIM(u) < DIM (v) 

RETURN u 
Loop 

If DIM(v) = DIM(u) exi t 
k_ :+ 1 

v := APPEND(v, [0]) 

Loop 

u := M0D(u - FIRST(u) ■ v , p) 
#22: u := REST(u) 

k_ :- 1 

If k_ < 0 exit 
v := DELETE (v , -1) 

Loop 

If FIRST(u) > 0 
RETURN u 
0 

RETURN u 
u := REST(u) 


Now I'm going to replace the routines polygcd(u,v,p,x) and polypower(u,k,v,p,x) above 
with the following new implementations. 
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po1ygcd(u, v, p, x, k_ := 0, t_, v_) := 

P rog 

u := polybovec(u , p, x) 
v := polytovec(v , p, x) 

If DIH(u) < DIM (v) 

[t_ := u, u := v, v := t_] 

Loop 

If V = [] 

RETURN vectopoly (u) 
t_ := INVERSE_MOD(FIRST(v) , p) 
v := M0D(t_- v , p) 

Loop 

If DIM (v) < DIH(u) 

P rog 
k_ :+ 1 

v := APPEND(v, [0]) 

#23: exit 

t_ := u 
Loop 

If FIRSTCO > 0 

t_ := MQD(t_ - FIRSTCO -v, p) 
t_ := RESTCO 
If k_ = 0 exit 
v := DELETE (v, -1) 
k_ 1 
u := v 
v := t_ 

Loop 

If V = [] 
exi t 

If FIRST(v) > 0 exit 
V := REST(v) 


3 2 

#24: polygcd(u, v, 7, x) = x + 2-x + 2-x + 4 

polypower(u, k , r, p, x , y , u_ := [1]) := 

P rog 

u := polytovec(u, p, x) 
r := polytovec(r, p, x) 

Loop 

#2S : If k = 0 

RETURN vectopoly (u_) 

If 0DD?(k) 

u_ vecrem(vecmult(u , u_, p) , r, p) 
u := vecrem(vecmult(u , u, p) , r, p) 
k := FL00R(k, 2) 

9S 76 54 3 2 

#26: polypower(u, 10, v, 7, x) = 2-x + 2-x + B-x + 4-x + B-x + x + 2-x + 6-x + 6-x + 1 

Again we got the correct result for our examples above and already now a noticable in- 
crease of performance can be seen, though the full power of these small gems of pro- 
gramming will become only visible when dealing with huge polynomials as the reader 
might again check himself. 

As for our goal, namely to provide all necessary tools for a later implementation of 
Schoof's algorithm, we need some further tools concerning computations on elliptic 
curves. The first routine of this sort computes the so-called division polynomials 
Y(n,a,b,p,x,y). They are of the form nf(x)y , if n is even and of the form f(x), if n is odd, 
where f(x) is some normed polynomial in x, whose zeros are exactly the 
x-coordinates of all n-torsion points (x,y) of E with y t 0 in the algebraic closure of Zp. 










D 


-N-L#67 


Surfaces from the Newspaper (2) 


p51 


Surface #2: ( x 2 + y 2 + z 2 ) 2 = x 2 + y 2 

Surface #2 in DERIVE presented applying functions from polycontour .mth (DNL#63) and 

IMPLICIT_Peter . mth (DNL#64). 

2 2 2 2 2 2 

VECT0R(ContourPts_2D((x + y + z ) - x - y , 1, -1, 1, -1, 1, 0.05, 0.05), 1, -0.5, 0.5, 0.1) 


2 2 2 2 2 2 


VECTOR (ContourDots_XY((x + y + z ) - x - y , 

1, -1. 

1,-1, 

LT, 

O 

LT, 

O 

1 

LT, 

o 
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2 2 

Imp! icitDots((x + y - 

2 2 2 2 
p z ) = x + y , [-1, -1, 

-0.5], 

[i, i, 

0.5], 0.05) 


2 2 2 2 2 2 
VECTOR (ContourPts_XY((x +y + z ) - x -y, 

1, -1, 

i,-i, 

1, 0.025, 0.025), ], -0.25, 

0.25, 0.05) 
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Surface #2 in DPGraph 



Surface #2 in MuPAD 

FA2 - Flache #2 

plot (plot : : Implicit3d ( (x A 2+y A 2+z A 2) A 2=x A 2+y A 2, 
x=-l. . 1 , y=-l. . 1, z=-l. . 1, 

Mesh = [21, 21, 21] , AdaptiveMesh = 3), 
Scaling = Constrained) 



Surface #2 in Autograph 3.20 



The Mathematics Education into the 21st Century Project 


together with 

The University of Applied Sciences (FH), Dresden (Germany) 

are proud to announce our 

10th (Anniversary!) International Conference 

“Models in Developing Mathematics Education ” 

September 11-17, 2009 
Dresden, Saxony, Germany 

in cooperation with 

Saxony Ministry of Education 


Chairman 

Dr. Alan Rogerson, International Coordinator of the Mathematics in Society Project 
(Poland). 

Prof Dr Fayez Mina, Dept, of Curriculum & Instruction, Faculty of Education, Ain 
Shams University (Egypt). 

You are invited to attend our project conference to be held in the historic and beauty- 
ful city of Dresden, Germany. 

The chairman of the Local Organising Committee will be Prof. Dr. Ludwig Paditz of 
the Dresden University of Applied Sciences. 

For ALL further conference details and updates please email arogerson@inetia.pl . 
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that will be announced shortly. 


